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A GENERAL THEORY OF THREE-DIMENSIONAL FLOW IN 
SUBSONIC AND SUPERSONIC TURBOMACHINES OF 
AXIAL-, RADIAL-, AND MIXED-FLOW TYPES 
By Chung-Hua Wu 


SUMMARY 

A general theory of steady three-dimensional flOT^ of a nonviscous 
fluid In siibsonic and. supersonic turhomachines having arhitrary hub and. 
casing shapes and a finite number of blades is presented. The solution 
of the three-d i mensional direct and inverse problem is obtained by 
investigating an appropriate combination of flows on relative stream 
surfaces whose Intersections with a z-plane either upstream of or some- 
where inside the blade rovr form a circular arc or a radial line. The 
eq^uations obtained to describe the fluid flow on these stream surfaces 
show clearly the several approximations involved in ordinary t\ro- 
dlmenslonal treatments. They also lead to a solution of the three- 
dimensional problem in a mathematically two-dimensional manner through 
iteration. The eqixation of continuity is combined with the eq_uatlon of 
motion in either the tangential or the oradial direction through the use 
of a stream function defined on the surface, and the resulting equation 
is chosen as the principal eq.uation for such flows. The character of 
this eq.uation depends on the relative magnitude of the local velocity of 
sound and a certain combination of velocity con^onents of the fluid. A 
general method to solve this eq^uation by both hand and high-speed 
digital machine confutations when the eq.uation is elliptic or hyperbolic 
is described. The theory is applicable to both Irrotational and rota- 
tional absolute flow at the inlet of the blade row and at both design 
and off -design operations. 


INTRODUCTION 

The problem of three-dimensional flow in turhomachines of axial-, 
radial-, and mixed-flos'j- types is treated in references 1 to 19. Because 
of the enormous mathematical difficulties Involved in the problem, Lorenz 
(reference l) first introduced the idea of an infinite number of blades 
of infinitesimal thickness in order to follow the flow on a given sintface. 
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Bauersfeld (reference 2) added to the theory the condition of integra- 
hility for the hlade surface that must he satisfied in the inverse, or 
design, problem. The theory is farther clarified and strengthened hy 
the vorks of Stodola (reference 3) , von Mises (reference 4) , and Dreyfus 
(reference 5), and is the basis of nsany recent investigations on axial-, 
radial-, and mixed-flow con 5 >ressors and turbines. 

For incon 5 )ressible flow, Ruden (reference 6) proves that the 
through-flow solution obtained under the assuniption of an infinite nim- 
ber of blades gives a circumferentially average value of the fluid prop- 
erties, provided the deviations of the fluid properties from their cir- 
cumferential averELges are small. In reference 7, Traupel points out the 
oscillatory nature of radial flmi in a multistage turbomachine and gives 
solutions of the three-dimensional potential flow through inclined 
stationary blades and also of the rotational flow through a homogeneous 
stage of Identical nontwisted blades for an incompressible fluid and an 
infinite numiber of blades bounded by cylindircal walls. Meyer gives a 
detailed treatment of three-dimensional potential flo\^ in a stationary 
blade row, for an incompressible fluid and cylindrical bounding -VTall, 
in reference 8, where the solution for an infinite number of blades is 
extended to a finite number of blades by the vortex^and-source method of 
Ackeret, which is originally given for two-dimensional flo\r (reference 9) . 
In reference 10, a linearized solution for an Incompressible fluid and an 
infinite numiber of blades for a prescribed loading and cylindrical 
bounding uralls is obtained by Marble, and is used later to investigate 
the problem of mutual Interference of adjacent bladfe rcn^s and off- 
design operations (reference 11) . Siestrunck and Fabri (reference 12) 
also obtained a linearized solution for incompressible flow, and the 
method is extended to compressible flow. For general wall shapes, 
Spannhake (reference 13) examiines the flow throu^ diffuser and impeller 
by the use of bound vortices for blades. The incompressible through 
flow in a mixed-flow impeller is treated by Gravolos (reference 14) . 

In reference 15, ¥islicenus examines the Influence on the meridional 
flow of the blade force and nonuniform circulation along the blade span. 

For compressible flo\7, Reissner (reference 16) gives a blade-design 
method in which the extension from an infinite number of blades to a 
finite mmnber of blades is accomplished by the use of a power series in 
the circumferential direction, and the terms in the series are deter- 
mined by a comparison of the eq^uations for an infinite number of blades 
and a finite numiber of bleides. (in reference 5, Dreyfus gives a method 
of designing water turbines of thin blades, in which the solution for 
an infinite mmiber of blades is extended to a finite numiber of blades by 
the use of a power series, the second term of which is determined from 
the equations of continuity and irrotational absolute flow and is 
explicitly given.) In reference 17 the compressible flow problems in 
axial turbomachines having an infinite number of blades are treated, and 
both the direct and inverse problems are considered. Methods for 
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limiting solutions for zero and infinite blade-row aspect ratios and a 
step-by-step method of solution^ as well as a simpler method based on an 
approximate knovrledge of the shape of the streamline, for a finite blade- 
row aspect ratio are given. Unaware of the_ work of Traupel at the time, 
the authors of reference 17 also en 5 )hasized the oscillatory nature of 
radial floi'r in multistage machines and suggested the use of a simple 
sinusoidal form of the streamline as a first approximate solution. Gdieir 
methods are derived for con^jressible flow, however, and are also extended 
to the case inhere both the hub ^d casing walls or either is tapered. 
Reference 18 gives a general through-flow theory for both direct and 
inverse problems and for subsonic or supersonic flo\f in turbomachines 
having arbitrary hub aind casing shapes. The supersonic through flow in 
rotating impellers having a prescribed flow along the casing and pre- 
scribed blade shapes is treated in reference 19. 

A general theory of three-dimensional flow in subsonic and super- 
sonic turbomachines of axial-, radial-, and mixed-flow types for a 
finite number of thick blades of finite thickness has been developed at 
the NACA Lewis laboratory and is presented herein. Both the direct and 
inverse problems are considered. The theory is applicable to either 
iorrotatlonal or rotational absolute flow at the inlet of a blade xavr 
and to both design and off -design operations. 

In the section BASIC AEROTHEEMODYMMIC RELATIONS, the motion and 
energy equations for the unsteady flow of a nonviscous compressible fluid 
in a rotating blade row are expressed in terms of the velocity components 
and of two basic thermodynamic properties of the fluid, namely, entropy 
and a modified total enthalpy for flaw in rotating blade rows with change 
in radial distance from the machine axis. Estimated entropy changes due 
to shock waves (in the case of supersonic flow), heat transfer (in the 
case of a cooled turbine) , or viscous effect can be easily accomodated 
in the calculation. The equations obtained show clearly the condition 
under which the flow through blade ro^fs can be treated on the basis of 
irrotational absolute flow. 

In the following section, a general potential equation is obtained 
for steady three-dimensional compressible flow through rotating or 
stationary blade rows when the absolute flow can be taken as irrota- 
tional. The methods of solution for both subsonic and supersonic flows 
are briefly discussed. 

A simpler method of solving the three-d±menBlonal irrotational 
(absolute) flow, which is also applicable to rotational (absolute) flow, 
is obtained by considering fluid flotrs on a numiber of relative stream 
surfaces whose intersection with a z-plane either upstream of or some- 
where in the blade row form a circular arc or a radial line. Equations 
governing the floi'r on these surfaces are obtained in the next four sec- 
tions. Q?hrough the use of a stream function defined on the stream 
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surface, the eq^uations of contlnialty and motion for fluid flow on these 
surfaces are combined into one principal equation. The character of the 
principal equation is dependent on the relative magnitude of the local 
velocity of sound and a certain combination of velocity components. 

The process involved in solving the direct and inverse problems hy 
this approach is described in the section STEPS FCER COMPLETE SOLUTIONS 
OF THREE -DIMENSIOML DIRECT AND INVERSE PROBLEMS. In the inverse prob- 
lem, besides the blade -thickness distribution determined by blade 
strength and other considerations, either the tangential velocity, a 
relation between the tangential and axial velocity, or one other rela- 
tion is prescribed on a mean stream surface about midway between two 
blades. The last section gives a general method of solution of the 
principal equation when it is elliptic or hyperbolic. 


SYMBOLS 


The foUovTing symbols are used in this report: 


a 


velocity of sound 


B,b 



C, c 


D_ 

Dt 




F,f 


H 

h 


integrating factor for continuity equation for S 2 and Sq 
surfaces, respectively 

differentiation coefficient used to multiply function value 
at point j to give the m"*^^ derivative at point i based 
on n^ degree polynomial 

nonzero term on right-haad side of continuity equation for 
S 2 and Sq surfaces, respectively 

specific heat of gas at constant pressure and volume, 
respectively 

differentiation with respect to time following relative 
motion of fluid particle 

m^^ derivative of q 

vectors having the unit of force per unit mass of fluid 

given function of on Sg 

total enthalpy per unit pass of fluid, h + ^V^ 

static enthalpy per unit mass of fluid, u + p/p 
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I modified total enthalpy for flow in rotating "blade row with 

change in radial distance from nachine axis, 

h + i CD^r^ or E - CD(Vur) 




k 

L 

M 


coefficients of first- and second-order derivatives in the 
principal eq.uation 

thermal conductivity 

distance along streamline 

orthogonal coordinates on surface of revolution 

mass flow het^reen mean stream surface and one surface of 
blade 


N number of blades 

n unit vector normal to relative stream surface S 


P 

Q 


1 


R 



r 



^2 


s 


T 


static pressure 

I 

heat added to fluid particle along its path per unit mass 
per unit time 

einy g_uantity on relative stream surface S 

gas constant 

remainder term of m^^ derivative at point i obtained by 
using n"^^ degree polynomial 

radius vector 

relative stream surface passing throu^ fluid particles 
lying on a circular arc upstream of or midway in blade 
row 

relative stream surface passing through fluid particles 
lying on a radial or curved line upstream of or midway in 
blade row 

entropy per unit mass 
s/r 

static tenperature 
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t time 

U velocity vector of tlade element at radius r 

u interval energy per unit mass 

V absolute velocity of fluid 

W velocity of fluid, relative to blade, V - U 

x,y independent variables 

z distance along turbomacbine axis 

^ l^u 
a arc tan — ^ 

r Wa 

X ratio of specific beats 

X average value of T for the teiiperature range Involved 

5 grid spacing 

6 equal to 1 and r for and. S 2 surfaces, respectively 

^ independent variable z or r for surface and z for 

Sg surface 

■q independent variable cp and r for and Sg surfaces, 

respectively 

9 angular distance of fluid particle measured -with respect to 

stationary radial line 

A slope of characteristic curves, U ^ 

X tan o 

. a 

|j, arc sin — 

"D equal to r and 1 for S]_ and Sg surfaces, respectively 

^ absolute vorticity, Vx V 

p fluid density 
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2 generalized varia'ble used for general density table 

a angle between tangent of streamline or boundary wall in the 

meridional plane and axial direction 

T radial^ axial, or angular thickness of stream sheet 

0 velocity potential 

0 generalized variable used for general density table 

cp angular distance of fluid particle measured with respect to 
radial line on rotating blade 

X angle between w and axial direction 

stream functions defined on relative stream surfaces S 2 and 
respectively 

00 angular velocity of blade 
Subscripts; 

c casing 

e exit 

h hub 

1 inlet 

I meridional component 

m mean stream surface 

0 lower limit of integration 

r,u, z radial, circumferential, and axial components 
s isentrbpic 

T total state 

T],^ components in t]- and ^-direction, respectively 

1 on S-|_, or in front of rotor 

2 on Sg^ 


or behind rotor 
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Supers crlpts : 

. . . k refer to points a,"b, . . . k, respectively 

BASIC AEROTHERMODYMMIC EELATIONS^ 

The three-dimensional flow of a nonviscous, compressible fluid 
through a turhomachine is governed, hy the following set of Basic laws 
of aerothermodynamics i From the principle of conservation of matter, 
the eq.uation of continuity is 


|| +V.(pW) = 0 


( 1 ) 


or 


V . W + 


D In p 
Dt 


(la) 


For a blade rotating at a constant angular velocity co 
z-eixis, Newton’s second law of motion gives 


about the 


+ 2a»xW *= - — vp 

Dt p 


( 2 ) 


Because the boundary walls are surfaces of revolution and the relative 
flow can be approximated as being steady in many cases, it is convenient 
to use a relative cylindrical coordinate system r, cp, and z with 
cp measured with respect to the rotating, blade (see fig. l) . By use of 


^-||+(W.V)W = ^.5 


-Wx(TXW) 


the scalar forms of the equation of motion (2) in the axial, clrcum 
ferentiai, and radial directions can be expressed as 

2 










dt 


u ^^u^u „ ^ 


u ^r^u _ 1 ^ 

+ + 2o¥^ 

r ^ pr ocp 




+ W- 


Wu ^z 






1 dp 

p ^ 


(2b) 

(2c) 


ISome of the relations given in this section have been given in refer- 
ence 18. They are repeated here for conpleteness and easy reference 
for the following developments. 
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The first law of thermodynamics may he written 


+ p d (p:3 .) = 


Dt 


Dt 


= Q 


where u is related to the tenqperature T hy 

du t= Oy dT 

and Q is given hy the following equation if only conduction is 
considered; 


(3) 


(4) 


Q = p"^ V-(kVT) 


(5) 


For the ranges of ten^erature and pressure encountered in ordinary 
turhomachines^ p, p, and T of the gas are accurately related hy the 
following equation of state: 


p = B pT 


( 6 ) 


Although the flow of the gas through the turhomachine is completely 
defined hy the preceding equations together with the Imown variations of 
Cy and k with tenperature and the given boundary and initial condi- 
tions^ it is found more convenient in references 17 and 18 to express 
the state of the gas in terms of the entropy and the total enthalpy or 
a quantity I of the gas, besides its velocity components. These 
quantities are defined as follows: 


T dfl « 
H = 

I = h + i 


du + p d(p"^) 

(7) 

b+|T2 

(8) 

- 1 u2 = H - co(V^;ir) 

(9) 


and 


h 6s u + pp“ 

From equations (lO), (4), and (6) is obtained 

dh = dT = Cp dT = ^ dT 


( 10 ) 

( 11 ) 


where y is equal to Cp/c^ and is a function of tenperature . Another 
expression for dh is obtained hy using equations (lO) and (7 ) , so that 

dh = ^ + T ds 


(Ua) 
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By the use of equations (7), (4), and (6), 




d In p - 


r-1 


d In p 


( 12 ) 


and 


r-1 


d In T - d In p 


(12a) 


can he obtained, and the equation of contintiity can he -written 

V.W+ 0 

^ r-1 Dt Dt R 


(13) 


Equation (13^ can he e^spressed in a slightly different form. From the 
definition of -fehe local -velocity of sound (reference 20) , 


(14) 


By the use of eq,uations (12) and (6), 


a^ S3 -f £ ta vET 

P 


(14a) 


Suhsti-fcuting this relation into equations (I2a) and (13), with the use of 
equation (U), results in 


- - dll . s 

dlnp = -^-dg. 

Qt 


(I2h) 


and 




(13a) 


From equations (9) and (Ua), 


-VP + ivW^ - CD^r S3 VI - Tvs 
P ^ 


The equation of motion (2) can then he -wri-fcten 

^ - WX(VX\HI) + 20)XW = - VI + Tvs (15) 

An alternative form of equation (15), which In-volves the vorticlty of 
the absolute motion, is obtained as follows: With the z-axis parallel 

to OJ, 
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hence 


But 


V = W + a»x r 


VX V =VXW +VX(a>X r) 


(16) 

(17) 


VX(o;Xr) = (r*v) o) -(wv) r + o)(v* r) - r(vo)) « 2w 


therefore 


VX V = VXW + 2ci) 


(17a) 


This relation can also he seen from the following expressions of rela- 
tive and absolute vorticlty expressed In terms of the rotating and 


stationary c ylin drical coordinates r 


(VXW), 


Cp^ 

1 

r 5cp 


/ V ^r 

(VXW)^ = ^ 


(vxw). 


1 

r 


V). 


(vxw) 


5r 

1 

'xW 


z and 
^u 

" -ST 

Stfz 

r <^cp 

■ sr 


u 


5F" " 3F 


, . 1 ^(Vu^) 


1 

r 


0. 


respectively: 


(18) 


J 


(19) 


and the relation 


3(v^r) a(¥^^r) 

Using equation (l7a) results in the alteniative form of equation (15) 


_ WX(VX V) = -VI + T7s 


(15a) 
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By use of eq^uations (2), (9), and (Ua), 


= + 

Dt Dt p Bt 


+ ^ (W*V) p + W • 


DW 

Dt 


U . 


Dt 


T 

Dt p 5t 


+ W-((D^r-2a»XW) 


- U 


DU _ Ds 1 ^ 

Dt " Dt + p ^ 


( 20 ) 


It may "be en 5 >haslzed that the preceding equation is a consequence 
of the equation of motion (2) and the thermodynamic relations (4), (6), 
(7), and (lO) . For steady relative flow, the rate of change of I 
along the streamline is seen to he proportional to the rate of change 
of entropy along the streamline. 

The energy equation (3) can he used to express the rate of change 
of entropy along the streamline hy the use of equation (7) as follows: 

S = ® i (21) 

The preceding equations lead to several iniportant general consider- 
ations: If the hlade rows are not placed too close together and no 

trailing vortices are shed from preceding hlade' rows (or where these 
effects can he neglected) , the fluid properties at a fixed point rela- 
tive to the hlade can he taken as constant with respect to time. Con- 
sequently, according to equations (20) and (21), the quantities s and 
I of the gas re main constant along the streamline for adiahatic flow. 
The invariancy of I means that the rate of change in total enthalpy 
along the streamline is equal to the angular speed of the hlade multi- 
plied hy the rate of change in angular momentum (about the machine axis) 
of the fluid particle along its streamline, which is the well-known 
Euler turbine equation usually derived under less general conditions. 

In a cooled turbine where the heat transfer may he large, the rate of 
change of s and I along the streamline can he corrected hy equa- 
tion (2l) for an estimated value of Q. Again, for steady relative flow, 
equation (l5a) shows that either when gradient I and gradient s both 
vanish or when the difference het:reen VI ^.nd TV s vanishes, the 
absolute vorticlty either vanishes or is parallel to the relative 
velocity. 

For the flow through a stationary hlade row cu = 0, W becomes V , 

I becomes H, and equation (l5a) becomes 

- V X (V XV) « - VH + TV s (15b) 

which agrees with similar relations previously obtained hy Vazsonyi 
(reference 2l) and Hicks, Guenther, and Wasserman (reference 22) . It 
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is interesting to see that, for relative flovr in a rotating "blade row, 
"VX(VXV) "becomes WX(vXV) and H "becomes I. 

If it is as sinned that the fluid enters the inlet guide vanes of a 
turhomachine with uniform H and s and zero vorticity and that the 
floif is adiabatic, s does not vary in the inlet guide vanes and p 
is then a function of only p, according to equation (l2). Conse- 
quently, by virtue of Kelvin’s circulation theorem, the a"bsolute vor- 
ticity will remain zero in passing through the inlet guide vanes and the 
flow in the inlet guide vanes can be treated on the basis of irrota- 
tional absolute flow. 

If the guide vanes impart a radial variation of tangential velocity 
of the fluid in a z -plane downstream of the vanes s i m il sir to that in a 
potential vortex, that is. Inversely proportional to the radius, the 
circulation is constant along the blade span and the fluid maintains a 
uniform s and H and a zero vorticity of absolute flow entering the 
following rotor -blade row. If the rotor -blade row is situated far away 
from the inlet guide vanes, the fluid enters the rotor with a uniform 
I in the circumferential direction, as well as in the radial direction, 
aufl the flow through the rotor blades can again be treated on the basis 
of zero absolute vorticity and steady relative flow. If the rotor is 
close to the guide vanes, however, vortices are shed from the inlet 
guide vanes because of periodic variation in circulation caused by 
unsteady flow, and the flow downstream of the stator and throu^ the 
rotor blades should theoretically be treated on the basis of rotational 
flow. 


If the guide vanes impeirt a radial variation of tangential velocity 
of the fluid at a z-plane downstream of the vanes not Inversely propor- 
tional to the reidlus, the circulation varies along the span of the guide 
vanes, vortices are shed from the trailing edge to the fluids downstream 
in the direction of the exit velocity, and the fluid enters the follow- 
ing rotor blades with a uniform s and E but a nonuniform I and a 
nonzero value of absolute vorticity. Consequently, the flow through the 
rotor -blade row can no longer be treated on the basis of zero absolute 
vorticity, even if it is far apart from the preceding guide vanes. 

From the preceding discussion, the choice of s and H or I 
as the two basic thermodynamic variables of the gas besides its veloc- 
ity conponents is apparent. Compressor and turbine rotors are usually 
designed to impart or subtract the same amount of energy to or from the 
gas radially; hence H is usually radially constant throughout the 
machine if the inlet flow is uniform (except in the boundary layer along 
hub and casing walls) . If the circumferential velocity of the gas 
upstream of the blade row is zero or varies Inversely with radius, I 
is then constant throughout the machine. These facts will be utilized 
in the following developments. 
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POTEHTIAL EQUAIION FOR THREE -DIMEKBIOKAL FLOW THROUGH 
ROTATING BLAKE RCW 


Consider first the special case of steady relative flow where the 
fluid upstream of the hlade row is free of vorticity and is uniform in 
I and s. The adiahatic flow through the "blade row is then relatively 
steady and absolutely imrotational and is most conveniently treated hy 
the use of a velocity potential $ based on the zero absolute vorticity 
and related to the relative velocity components throu^ eg.uation (16) as 
follows: 






1 

r ^ 


"Vu =» Wu + on* 








> 


J 


(22) 


For steady isentropic flow, the continuity eq^uation (l3a) becomes 


1 ^(W^r) 
r 5r 


1 


dz 




Sh 

SF 


W-- -V- + 


^u Sh 
+ — ^ + 




Sh 


:)=o 


(23) 


From eq.uations (9) and (22), 


T ^ 


3h 


^2$ 




1 Sh 
r ^ 



dh 




,Wu ^25, Vu2 

r dr^ ” r ^ 2 ; dzdiy 


drdcp ^ ^2 ^2 ^ r dcpdzy 


dzdr r dcpdz 


+ ¥, 



(24) 

(25a) 


(25b) 


(25c) 


By the iise of equations (22) and (25a) to (25c), the continuity equa- 
tion ( 23 ) may be written 



3V 


MCA OT 2604 


15 


/ 32 $ / W„ 2 \ 32 $ / w/\ 32 ,, 

^ ■ aV 31.2 V aVrW \ ■ ^ 2 / 3^2 - 

WjWq 32 $ ^uVz 1 32 $ 32 $ j. 1 ^ j. Z?L_^ n 

0^2 r (^rdcp “ a^ ^ dT6z " g^2 3z5r r y ^ J "Sr 

(26) 

Equation (26) Is then the three-dimensional potential equation for 
Isentropic flow in a rotating blade row. It is seen from this equation 
and equation (16) that the three-dimensional flow through a rotating 
blade row cannot be treated by a three-dimensional flow throu^ a sim- 
ilar stationary blade row with the same inlet condition relative to the 
blade row, as in the case of two-dimensional flow on a cylindrical sur- 
face, because the difference between the absolute and relative vortlcity 
2oJ does not enter into the two-dimensional flow on a cylindrical surface 
but does enter into the three-dimensional case. 

Equation (26) is veiy similar to the ordinary three-dimensional 
potential equation for flow past stationary objects, except that both 
relative and absolute velocity coniponents are involved in the coeffi- 
cients of $ derivatives and that is directly defined by the abso- 
lute velocity. The real difficulty in solving this equation lies in 
the fact that all the velocity congjonents change greatly in passing' 
through a turbomachine and, consequently, the equation cannot be 3Jn- 
earized and yet give a good approximate answer. For supersonic relative 
velocity, the method of characteristic surfaces (references 23 and 24) 
maybe used to solve equation (26), with the initial conditions not 
given on a characteristic surface. For subsonic relative flow, the 
equation is more conveniently written in the form 


1 3$ 1 3 In p ^u S In p 


+ 


0 

oz 


(26a) 


and can be solved by Southwell's relaxation method (reference 25) or 
other numerical methods using the differentiation formulas obtained in 
reference 26 to take care of the unequal grid spacings near the blade 
surfaces and the curved hub and casing walls. The last three terms in 
equation (26) are computed from the values or velocities obtained 
in the previous cycle and kept as constants during each Improvement of 
$ values, and the whole process is repeated until the desired accuracy 
is obtained. Because a three-dimensional stream function cannot be 
defined, the use of velocity potential results in a boundary-value prob- 
lem of the second kind, which is more difficult to handle in the calcu- 
lation than the first kind. The bcwndary condition to be satisfied is 
that the relative velocity normal to the miovlng blade is zero, or 
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Wj. rij, + Wjj +’^2 ® 


(27) 


where n is the unit normal vector at the houndary surface, and that, 
at inlet and exit stations far away from the hlade, the velocity is 
parallel to the hounding huh and casing walls, which, in the case of 
the axial machine, means that 


Vr = ||=0 at 


±c» 


(27a) 


In hoth the subsonic and supersonic cases, the solution is 
extremely time-consimiing. Furthermore, this direct approach to the 
three-dimensional problem requires that the absolute velocity at the 
inlet to the blade row be irrotational and the flow be adiabatic. In 
actual mchines, the flow entering the blade rows is always rota- 
tional, -vdiich is caused by a nonuniform total enthalpy and entropy at 
the inlet of the nachine, by entropy change caused by shock waves or 
heat transfer, or by the effect of boundary layers along the hub and 
casing walls. Some other approach to the problem, which is simpler to 
handle and is also applicable to rotational inlet flow, is therefore 
desirable. One approach is suggested in the following sections. 


FOLLOWING FLUID FLOW ON RELATIVE STREAM SURFACES 

In order to solve the steady three -dimensional flow, with either an 
irrotational or rotational absolute flow at the inlet, in a relatively 
simple manner, an approach is taken to obtain the three-dimensional 
solution by an appropriate combination of mathematically two-dimensional 
flows on essentially two different kinds of relative stream surface 
(figs. 1 to 3) . The first kind of relative stream surface is one whose 
intersection with a z-plane either ipstream of the blade row or midway 
in the blade row forms a circular arc (fig. l) . The second kind of 
relative stream surface is one whose intersection with a z-plane either 
upstream of the blade row or somewhere inside the blade rov?- forms a 
radial line (fig. 2). These two kinds of relative stream surface will 
be hereinafter designated stream surfaces Sq and S2^ respectively. 


S]_ Stream Surface of First Kind 

In figure 1 is shown a stream surface of the first kin d formed by 
fluid particles lying on a circular arc ab of radius oa upstream 
of the blade row. It is usually assumed in ordinary two-dimensional 
treatments (for example, references 27 to 30) that the stream surface 
thus formed is a surface of revolution. In the following development, 
the surface will be allowed to take whatever shape it should have in 
order to satisfy all the equations governing the three-dimensional flow. 
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In most cases, the deviation of the surface from a surface of revolution 
is not large, and it is satisfactory to consider S-j_ surfaces formed hy 
fluid particles originally lying on a circular arc upstream of the blade 
row. If the rotatlonality of the inlet absolute flow is large, if the 
blade is designed for a velocity diagram quite (different from the free- 
vortex type, or if the blade length is long in the direction of the 
through flaw (radial- and mixed-flow machines), the twist of the surface 
may be quite large, resulting in very large circumferential derivatives. 
If this effect is found during calculation or known from experience, it 
is more satisfactory to consider S 3 _ surfaces formed by fluid particles 
originally lying, in front of the blade row, on carves inclined to the 
circular arc in a direction opposite to the twist of the surface. In 
this way, the Intersection of the S 3 _ surface with a constant z-plane 
about midway in the flow path is a nearly circular arc, and the total 
twist of the surface will be about equally distributed toward the 
iqjstream and downstream directions (fig. 3) . If this distribution of 
the twist of the stream surface is still not enough, it may be necessary 
to divide the complete flow path into a few shorter paths and consider 
an S]_ surface for each of them. Under these conditions, Sq surfaces 
formed by fluid particles originally lying on the hub or casing walls 
upstream of the blade row should not be chosen in order that the congoli- 
cation arising from the possibility of fluid particles leaving the wall 
and flowing along the blade surface may be avoided. In such cases it 
is better to consider the Sq surface a short distance from the hub^^and 
casing; otherwise, for an approximate solution the fluid can be con- 
sidered to follow the hub and casing walls, which are surfaces of rev- 
olution, and the calculation is thus much simpler than that for a 
general surface. 


S 2 Stream Surface of Second Kind 

A stream surface of the second kind is shown in figure 2. The 
most important surface of this family is the one about midway between 
two blades dividing the mass flow in the channel into two approxinstely 
equal parts. This surface is designated the mean stream surface 
(Sg,m) • blades id.th radial elements, such as the one shown in fig- 

ure 2, it is convenient to consider a mean stream surface formed by 
fluid particles originally lying on a radial line ab upstream of the 
blade row if the twist of the surface is not expected to be large. 
Otheiwise, the radial line is chosen about midway in the passage with 
the fluid particles originally starting out from a curved line iqistream 
of. the blade row such as shown in figure 3. 

The mean stream surfaces for axial-flow gas turbines designed on a 
free -vortex velocity diagram are shown in figures 3 and 4. The radial 
element of the mean stream surface (fig. 4) is, chosen accordingly as the 
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stator is designed to aline tte "blade sections radially at the leading 
edge, trailing edge, or somewhere between. Inasmuch as the rotor -blade 
sections are. usually alined radially at or near the center of gravity 
of the blade sections, the radial position of the mean relative stream 
surface is chosen at the same position (figs. 3 and 4). The continua- 
tion of the stream surface outside the blade row is not shown. The mean 
stream surfaces for the inlet stage of a multistage axial compressor 
designed on the principle of a symmetrical velocity diagram at all radii 
are shown in figure 5. 

Both of these two kinds of stream surface are employed, in general, 
in the solution of the three-dimensional problem. The correct solution 
of one surface often requires some data obtainable from the other, and, 
consequently, successive solutions between these two are involved. Yet, 
the solution of each surface is manageable with the present mathematical 
technique aind conputational facilities. In many practical cases, and 
especially in the inverse problem, however, this iteration may not be 
required if only an approximate solution is required or if the prescribed 
values lead to a satisfactory blade shape. These points will be dis- 
cussed in the section next to the last (pp. 53 to 57), 


Relations among Relative Velocity of Fluid, Coordinates of 
Stream Surface, and Notmal to Stream Surface 


In general, the coordinates of the stream surfaces and their differ- 
entials axe related, respectively, by the following equations; 


S(r,cp,z) » 0 


as 




as 


(28) 

(29) 


Rather than use the three partial derivatives of S with respect to 
the coordinates, it is convenient to consider the unit vector n normal 
to the surface, which is related to S by 


°r 

a? 


as 

r a^ az 



(30) 


The vector n is, of course, perpendicular to the relative velocity W, 
so that 


n .W =» 0 



MCA IN 2604 


19 


or 


+ V^u + = 0 (31) 

By using equation (30 ) , eq.uation (29) can "be -written 

nj.dr + n^rdcp+n 2dz = 0 (29a) 

The vectors n and W are shown on S]_ and 83 surfaces in fig- 
ures 1 and 2. 


EQUATIONS GOVERNING FLUID FLOW ON SURFACE FOR AXIAL-FLOW AND 

AXIAL-DISCHARGE mXED-FLOW TURBOMACHINES 

If the fluid motion on S^ is followed, equations (28) and (31) 
can he used to e1.imi.nate one of the three coordinates. For axial-flow 
(figs. 6(a) to 6 (c)) and axial-discharge mixed-flow turhomachlnes 
(fig. 6 (d)), it is convenient to express r in terms of Cp and z. 

For radial -flow and radial-discharge mixed-flow turhona, chines (figs. 6(e) 
and 6 (f)), this system -will encounter difficulty at the exit where the 
rate of change of fluid state -with respect to z becomes infinite. It 
is therefore necessary to eliminate z and to consider r and cp as 
the two independent -variables. 


Flow Along General Sj^ Surface 

For axial- and mixed-flow turbomachines, any quantity q on the 
Sj. siarface is considered a function of cp and Z} that is, 

q = q[cp,z,r(cp,z)] 


The change in q along Sj_ due to a small dhange in qj while z is 
held constant is (see fig. l) 


dq t= 



From equation (29a), 
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hence 




°r 


Slmilaxly, for dcp = 0, 


dq. 


'Sq. _ % aq\ 

,^5z ~ 5r/ 


dz 


With a hold partial derivative sign used to denote the rate of 
change of any quantity q. on "with respect to q> or z, with the 

other kept constant, the preceding relations give 


.1 ^ 1 Sq. ^ Sg 

r 3cp r “ Uy Sr 


la 

3z 




(32) 


With the ^relations (31) and (32) , the rate of change of cl along a 
streamline on'' St is 


2a _ 2a + u la 

Dt " r aqp "z 9z 


(33) 


Eq^uations of continuity and motion. - When the fluid motion is 
followed along the stream surface and equations (31) and (32) are used, 
the continuity equation for steady relative motion becomes 


where 


c(<p,z) 


1 

r 


9(pWu) 3(pwj 

3cp 3z 


1 (°r 
Uj, \r ^ 


= P c(cp,z) 

+ % + “z 



(34) 

(35) 


For rotational steady relative motion, the equations of motion (l4) 
in the radial, circumferential, and axial directions are 
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Relations (9), (l 6 )^ (3l), and (32) along the relative stream sur- 
face Si can be used to reduce equations (36) to the following: 


Wu^ Wu 9Wr aWr 

r r acp 3 z 


rp ^ , ' 2 ^ 


A 




- W-r 


aWj 

az" 


/l ^^z _ 1 31 T as ^ (hh _ Ss 2 \ 

\r acp az/ ■*■ - "r acp+r 

. /l ai „ as ^Z /hh _ o \ 

^r w ■ az y “ • az ^ ^ ^ ^ 


(37) 

The last term in each of the preceding three equations is propor- 
tional to the conponents of the normal vector and therefore can be 
expressed an a component of a vector that is parallel to n and has 
the di m en s ion of force per unit mass. If this term is defined as 


f = 


_ /Sh 

n^ \3r 


T 


^s 



n 



(38) 


the preceding equations can be written 


^u 

r acp ^z az 



- 2 oW^ = f 3 , 


(39a) 


WrWu 


Wj, aWj, 




1 '’"z 

r acp 


+ 2 clWt 




aw^ 

az 


+ w, 


u 


3Wu\ 

8z J 

/ 1 ^ 

\r acp “ az / 


1 OT T ^ 

r acp r acp ^u 


lro'+ fu (39b) 


+ T 1^ 
9z 3z 


+ fry 


A 


(39c) 
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Similarly, the equations of motion, in. the form of equations (2) can he 
■written 


it I TJ ^ j=2 -P * 1 

r acp ^ 3z r 


Wu. aWu. 
“ 8^ 


+ Wz 


3Wu. Wj-Wu. 
3 z ^ "r 


+ 


2a^Ir 


— — + f ' 
pr aqp u 


^3¥^ 
r W 


+ W„ 


aWz 

3z 


1 3p 

p 3z 


+ f ’ 




(40) 


•where 


f = 


-L i ^ 

n^ p ^ 


Because this vector f is parallel to n , it is peipiendicular to 
the relative velocity of fluid, or 


frWr + fuWu + f z^z = 0 


(41) 


By the use of equation (41) and equations (39), it can he shown -that for 
steady flow on an Sq_ surface. 


DI „ Ds 
Dt Dt 

Which agiees -with equation (20) . Therefore, for the present problem of 
steady relative flow on a s-fcream surface, the relation (41a) can he 
taken either as one of -the equations of motion or to represent the rela- 
tion given hy equation (41) . In other -words, there are only four inde- 
pendent relations among equations (39a), (39h), (39c), (41), and (41a). 

Just as in the case of the continuity equation, ei-ther set of the 
preceding equations of motion is e3q>ressed in tenns -of. the special par- 
tial deri-vatives -with respect to the 'two independent -variahles qj and 
z. The effect of radial pressure gradient is -taken into account in all 
these equations hy the f term, which is neglected in the ordinary -two- 
dimensional -treatment on a surface of revolution. Equations (28), (31), 
(34), and (39) or (40), however, lead to a possihillty of correctly 
sol-ving the three-dimensional flow of fluid particles on an S 3 _ surface 
in a mathematically -two-dimensional manner. 

Principal equation. - The equations of continuity and the equation 
of motion in the circumferential direction can he combined into a prin- 
cipal equation through the use of a s-feream function ijr as follows; 
First, if a -variahle h is introduced such that 
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or 


D In b . 

Dt ~ °r ^ 


(42) 



^z ^z\ 

“r r/ 


dx 



(42a) 


in vMcb the integration is performed along a streamline on tbe S]_ 
surface, then the continuity equation (34) can be -written 


a(bp¥ ) a(bpW r) 

+ t3 0 


aqp 


az 


(34a) 


The preceding equation is the necessary and sufficient condition "that 
there exist a function ilr -with 


II = rbpW^ (43a) 

i = -1>PW„ (43b) 


The difference in ilr^ at -fcwo points J and k on the surface is 



•'k 

bp(W2.r dcp - 


\J 






In particular, the difference in ijr at -fcwo points j and k on the 
constant-z plane at the inlet -where the fluid s-fcate is uniform is 





These two equations show that, physically, the Integrating factor b 
can be interpreted as proportional to the local radial thickness of a 
thin stream sheet whose mean surface is the s-fcream surface considered 
here. The continuity equation (34a) can also be ob-balned by consider- 
ing the mass flow going into an element of such a s-tn*eam sheet as shown 
in figure 7. By equating -fco zero the mass flow going into the element, 
which is defined by -fcwo axial planes dcp apart and -fc\TO normal planes 
dz apart (see fig. 7(a)), and letting dcp and dz approach zero, 
there is obtained 
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3 (tp¥^) aCrpw^r) 
az 


(34b) 


where t Is the radial thiokness of the stream sheet. From equa- 
tions (34a) and (34b ) , it is apparent that b is proportional to T, 
and the differences in ijr at two points j and k as given by the 
two eqimtions preceding equation (34b) are proportional to the mass flow 
across any line joining the two points. Jn actual computation, only the 
ratio b to b^ or T to Tj_ is inportant (a different initial value 
amounts to a different constant multiplier of the relation between ijr 
and mass flow) . In the following, b will be retained in the equation, 
but in actual calculation it is simpler to evaluate the ratio T to 
than to evaluate the ratio b to b^, both from the data obtained 

on the Sg surface to be discussed later. Although the evaluation of 

this ratio requires, in general, calculations on the Sg surfaces, a 
means is nevertheless provided to deteimlne correctly the flow on a 
general S-^ surface throu^ itearation. 

From equation (43), 


■h 1 _ 1 a^ijf 1 a In bp at . “u i aijr 

- ~ " o - ^g gcp acp”** 3,2 acp 


r acp 


r“ acp^ 


bp 


a¥, 


u ^ a In bp a-ijf 


az 


3Z 


az 


(44a) 


(44b) 


The third terms in the preceding equations can be expressed in terms 
of h through the use of equation (l2b) : 


dlnbp=sdlnb + — g- dh - ds 

a 

where s* = s/R. But from equations (9) and (43), 

2 r, ,2 


(45) 


h I + 


OD 


,2^2 


2 2 

Then from equations (45) and (46), 


-(S)J 


\j? - (^fu 




)] i - i A vn- 


03^r^ " 


acp 


r acp 


a In b 
acp 



(■bp)' 




-•2 ^ a^paz^ 


°r 
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Substituting the preceding two eq^uations into equations (44a) R.nfl (44b) 
and adding yield 





+ 




1 3 (l , ^ 3 In t 33*^ 

_r acp \ 2 / r \ acp “ acp/ 


a2-w^2 

r • n^ 


[1 

r acp 



S^stituting the preceding equation into equation (39b) and dividing by 
a give the principal equation for the determination of fluid motion 
along a general S^^ surface; 


1 


W. 


u 


r^ aq? 


?!+_ + ?!i , N 3+ 3i|r 


where 


(48) 


az ^ az ^ g2 


31 aWj, cD^r^ + n„ 

az ^ az r n^ 


1 a In b 1 as* 1 ( 1 al Wr aWr a2_(Wu^-Wz^) Wr 
r acp "^r acp "^^2 \ r acp’*' r a^'*' ? W“ + 

' Z 



i §I +?. 9s 
r acp r acp ^U + 




(i 

\r acp 



The equation of the characteristics of the differential equation (48) is 
(reference 31) 
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+ 


2 — ^ r 


^4. 

dz 



from vMcli 


r 


dz 


a2-w/ 




(49) 


(50) 


Eq.uation (50) shows that the characteristics are real vhen 

+ ¥ 2 ^ > a, in which case the method of characteristics for two inde- 


pendent variahles (references 20, 30, 31, and 32) can he applied. When 


/\|^ 


+ Wz^< 


the characteristics are imaginary, and it is more con- 


venient to solve the eq^uation hy releixatlon (references 25, 33, 26, and 
29) and matrix methods ("references 26 and 29) in the following form, 
which is obtained hy suhstltuting equations (44) into equation (39h) : 


1 a^ijf a^T{f 1 / a In hp °^u\ aijr a in hp 3i}f 

3cp^ 3z^ _r^ \ 



1 31 ^ T as ^ ^ ■^^r 

racp racp u. racp"r 



0 (48a) 


Procedure of solution. - It may he noted that equation (39h), 
instead of (39c), is chosen to form the principal equation (48) or (48a), 
because f^j is, in general, much smal.ler than f 2 » The various quan- 
tities appearing in equation (48) or (48a) are to he computed from 
other equations given earlier o With the introduction of the stream 
function, there are adtogether seven basic independent relations - one 
energy equation (21) j three equations of motion, (48) or (48a), (39a), 
and (39c) j two equations between ijr derivatives and fluid properties, 
(43a) and (43b) j and the orthogonal relation between W and f, equa- 
tion ( 41 ) or (41a) . On the other hand, there are ten basic dependent 
variables in ijr, b, W^., W^^, ^r-’ ^u^ ^z^ and I (or h or p) 

to define the flow and the shape of the surface. In general, the vari- 
able b is to be evaluated according to equation (42a) or from the var- 
iation in the radial thiciaiess of stream sheet using the data obtained 
in the solution of S 2 surfaces euid is therefore considered as given 
here. If during the complete solution of the three-dimensional flow the 
shape of an S]^ surface is taken as the one obtained by Joining corre- 
sponding streamlines obtained on 83 surfaces of the preceding cycle, 
two relations between the n- (or f-) components are given by equa- 
tions (29) , and there are now altogether nine equations to be solved to 
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find the nine unknowns . AltermtlTely, the variation of may he 

considered as known ffom the S 2 solutions of the previous cycle^ and 
the remaining eight variables ^ which determine the flow on and the shape 
of the S]_ surfaces, can be determined from the seven preceding rela- 
tions given and the following additional relation; Because f^, f^^, 

and fg, respectively, are proportional to and ^ of the 

integral surface S, they satisfy the following eq[uation (reference 34) 

f .V X f « 0 (51) 

which may be written 


1 


1 

r 5r~l 



+ f-: 


1 

r S* 


1 

r dcp_ 


0 


(51a) 


By using equations (3l) and (32), equation (51a) becomes simply 

A - _8. 

9 cp az Vfr / 


(53b) 


This equation can be used to give fy_ by integrating along a constant 
<p line; 



(51c) 


If at z S3 z, 


o’ 


f = 0, then 


'•u 


•"r 

r 


acpVf, 


dx 


(51d) 


3In this case, then, the shape of the S]_ surface is determined after 
the f-conponents (or n-con 5 )onents) are obtained in the solution. In 
either case, equations (21) and (41a) are Invariably to be used first to 
determine the change of s and I. If the flow is isentropic, s 
and I remain constant along its streamlines on the surface. (For 
such a case and for a uniform Iriet condition, p in the continuity 

Y *— 1 

equation may be replaced by h and, consequently, the fp and z 
derivatives of s, as well as those of I, will not be involved in 
the equations (46) to (48)). In case of heat transfer or shock, the 
changes in s and I can be estimated by whatever method is available 
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and used In the calculation. For supers ord.c flow^ all the eq.uations are 
used to coD 5 )ute the fluid state at each point and the solution is carried 
downstream step hy step. For suhsonic flow. Iteration over the whole 
domain is necessary. The details of these computations will be given in 
the last section. In general, the solution of the flow on the general 
S]_ surface is very laborious, and is to be used in the final stages of 
calculation of the complete three-dimensional problem or when a high- 
speed computing machine is available. 

If the flow is such that it may be assumed to take place on a sur- 
face of revolution (at the hub and casing walls or other radii), the 
equations are considerably simplified as follows; 


Flow along Siarface of Revolution 


When the S]_ surface is a surface of revolution. 


, n^ = f ^ = 0 


(52) 


Let 




Wr 

=s tan u = A 


(53) 


where A is a given function of z. (For a conical flo\'?’ surface, A 
is simply a constant.) Equation (35) now gives 



(54) 


^•Thether c can be taken as zero will be determined by the relative 
magnitude of the three terms on the right side of the equation. In 
general, for nonnegligible c, equations (43) now become 

^ = r'bpWz (55a) 
II = - bpW^ (55b) 


Because Wj. is now related to Wz by equation (53), the three velocity 
components . can be solved simultaneously as follows : By use of the rela- 
tions (52) and (53), equation (39b) can be changed to 


1 X . p'] 1 /l 

(1+^ ) ? ^ - w; Vr 


hi 


T Ss 
r. 39, 


= 0 


(56) 
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Instead of equation (46) , 


h = I + 


aj2r2 




should he -written. Then 


+ («xW) g - 


(t»p) 


2s a^Tjf 
r 6cp3z 




- ¥. 


a^\lr 

%z2j 


(58) 


Combining the preceding equations -with equations (55) and substituting 
the resulting equation into equation (56) give the following principal 
equation for -fche flow on a surface of resolution; 


(1+A^) 


) -i. 


(1+X^) 


WuWz zh (. ^r^+^z^^ 

a2r a2 / 


¥ Silr 
az2 r ^ 


+M 


3Tlf 

dz 


= 0 


(59) 


where 



9 In h ^ 8^ 

9z 3z g^2 



a2 -w^ -W_2 -w 2 _o)2r 2 

X Zi 

r 



N = - (l+A^) 


1 a In b 1 9s* 1 1 Si 

a2-¥2 

2coX 1 /l Si T Ss\ 

r dcp r Sep g^2 r 

a2 

^z ^¥/ \r S^“r ^ 


For this equation^ the characteristics are real or imaginary when the 
resultant relative velocity ¥ is supersonic or subsonic, respectively. 
For the subsonic case, it is again better -to use the follo^rf.ng form 
ob-fcained by differentiating equations (55) and substituting the result- 
ing equations into equation (56) ; 
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- 2>(xibp - 


(bp)^ ( 1 Si _ T ^s\ 
1 d\lf \r ^ r 


(59a) 


With X given and h determined from data obtained on the S 2 
surface, there are now the six independent relations eq^uations (21), 
( 59 ) or (59a), (55a), (55b), (53), and (41a) for the determination of 
the six main variables in ijr, W^^, \t^, W^,, s, and I. 0?he f-components 

are not involved in the calculation. If the flow is adiabatic with 
uniform I and s, the equations are further sin 5 >lified. 


Flow along Cylindrical Surface 

If the flow near the walls of an axial-flow turbomachine can be 
considered to take place on a cylindrical surface, then 

^u. “ ~ ^ ^z ~ ^r “ ® (so) 


Equation (35) now gives 


c » — 


SWr 

SF~ 


(61) 


which is relatively small, (if c is negligible, b can be taken as 
1 everywhere.) For flow without change in radial distance, the quantity 

I can be replaced by (^ = h + The equations governing the 

cylindrical flow are then (conq>are reference 29) 


% = 

(62a) 


(62b) 


(63) 

® ^ „ Ds 

(64) 

171 ^ 

Dt " Dt 
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and 


1 - 


¥ 


u 


1 


^ _ 2 


j 5cp^ 


1 '\ ¥ ^i}r dijr 

^ (65) 


vhere 


,, ^ In t) 3s* 1 

M= ^ + 2 ^ 

a 


dz ' 3z 


1 3 In "b , 1 3s* , 1 
“ " X 3^^ r 3cp g^2 


a^-Wu^ 1 3h^ a^-w^ T 3s ' 

2 r ^ ^^2 r ^ 


W, 


or 


A. 


/ 1 3 In bp 3 

4 3 In bp 

3\}f\ 

("bp)^ / 


T 3s\ 

r^ 3cf? 


i^j,2 3 

5z 

33y 

\ 

r 3cp 

■ r 3^y 


(65a) 


In general, the circumferential derivatives of and s are to 

he determined hy the inlet flow and equations (63) and (64).. For adia- 
batic flow with uniform Hy and s upstream of the blade row, these 
derivatives are equal to zero everywhere, making the problem much sim- 
pler, The main dlfferehce between this simplified case and the ordinary 
t\ro -dimensional flow on a cylindrical surface is the inclusion of the 
factor b in equations (62) and (65) (in general, b is a function of 
cp and z) . If the velocity diagram is such that there is considerable 
radial gradient in the radial velocity or considerable variation of the 
dlstajice between the adjacent streamlines, the factor b is not 
negligible. 


EQUATIOMB GOVERNING FLUID FL0¥ ON S]_ SURFACE FOR RADIAL-FLOW 
AND RADIAL-DISCHARGE MIXED-FLOW TURBOMACHINES 
Flow along General Sq Surface 

For turbomachines with radial discharge, r and cp are considered 
as the two independent variables j that is, 


q = qj^r, cp, z(r,cpj] 
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tpiien 


3q. _ ^ ^ ^q. 

3r ~ 5r “ ^ 

1 ^ S3 2: ^ ^*3- 

r acp r ~ rig ^ 

5a - w M 5a 

Dt "1“ 3r ^ r acp 




J 


( 66 ) 


Equations of contlntilty and motion. - By the use of these relations, 
the equations of continuity and motion become 


a(pWj^) a(pw^) 


ar 


+ — 
r 


•where 


kk 




acp 
SiW, 


pc ' 




^"Z] 


(67) 


(68) 


and 


W. 


.2 

u 

r 


- W. 


Wr^u 


+ W 



aw^, 


ar 

r acp 

" ^Z ar 

aw^ 

Wj, aWy 


■ ar 

■ IT acp 

r aqj 


- 2oW^ = - g + T I? + (69a) 


i 55 4- 5 5a + f I 

r acp r acp u 


•with 


DWz aw w^aw^ 

Dt ar r acp 


1 /dh Ss\ 1 1 Sp 

n^ \3z “ "^z/ " “ n^ p 5z " 


(69b) 

(69c) 

(70) 


Principal eq-uatlon . - If a ■variable b' is introduced such that 

D In b* 


Dt 


= - c ' 


(71) 
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or 


1 a 


PL 


C* dX ra 




(71a) 


in which the integration is performed along a stream lin e on the surface, 
then the continuity eq.uation (67) can be written 


a(b«pW^r) 3(b«p¥^^) 


0 


ar acp 

and a stream function > can be defined on the surface with 


(72) 


3r 


-b*pW^ 


(73a) 


rb-pWr (73b) 

Here b* can be interpreted as .the thickness (in the z-direction) of 
the stream sheet whose mean surface is the S]_ surface considered. The 
continuity equation (72) can again be obtained by equating the mass flow 
into €ind out of an element of the stream sheet as defined by two axial 
planes dcp angle apaxt and two cylindrical surfaces dr distance apart 
as shown in figure 7(b). As before, the difference of ijr at any two 
points on the Sq surface is equal to the mass flow across any line 
connecting these two points. By the use of the preceding two equations 
and the relation 



the principal equation for the flow of this surface is obtained from 
equation (39b) ; 



WrWu 1 a2ijf / 

g2 r 3r3Cp y 



_1_ 

r2 aq? 


+ M 


3r 


+ 


N 31^ 
r ^ 


(75) 


where 



3 In b * 3s * 1 

3r 3r q2 


3^2 a2-(¥/l¥u^) cD^r^-Wp^ 

+ ¥2 ^ + 

3r ^ 3r r r 
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W = - i 


h» . 

1 as* 

+ 1 (_ 

91 + 

acp ■*' 

r acp 

a2r \ 

acp 


f_ + 

T as 

aV 

2 

\ r acp + 

r aqj 


aw. 


+ f ’ 
^ u 


+ 


20^3, + 


Wz 9Wz\ 
r acp/ 


This equation is seen to he hyperbolic or elliptic when 

greater or less than the speed of soxind, respectively. For the elliptic 

case, it is preferable to use the following form: 


a^ a^Tjr ( a in h*p a-ijf a in h 'p aV\ , 

a^2 3^2 a^p2 " V ar ar 3,2 acp aq^/ 


1 a^ 
r acp 


r acp r acp + r acp 




= 0 


(75a) 


The integrahility condition (51) is now written 


hence 



(76) 



(76a) 


The procedure of solving the principal equation ^■lith the various terms 
in it determined by other f lovr equations is the same as that in the 
previous system. 


Flow along a Surface of Revolution 

For the special case of flow on a surface of revolution, equa- 
tions ( 52 ) and ( 53 ) hold (with ?\ considered as a function of r) and 
the expression of c’ reduces to 


( 77 ) 
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Furthermore, equations (73) hecome 


-k'pWu 

rh'pW^ 


and equation (69h) becomes 


1 + 


I SWj, aw. 

■\2j r 3r 

Using the relation 




^ ^ I 1 fl^l T Ss\ 

- w^ = 0 


h = I + 


cu^r2 




gives the corresponding principal equation as 




) '’/l I 1 A I ^ 1 


vhere 


M = 


a In b« as* 1 

ar ar gZ 


ai ^ a2-w^-w^2-w/-<u2r2 w^2 

9 r r " ?v dr/ 


W = 
or 


. 1 Vl ^ In b' 1 as* . 1 ai\ a2-w2 

\ >?h~^ a^r^/ 



Wr^r 



as 


a^Tlf . 1 ^ l\ 1 fs In b'p a\lr 

ar2 + r ar + + ^2j ^2 "[_ ar ar 



1 a In b'p 

^2 3^ 


2 oib 'p 




'i ai _ T 

3 ^ r acpy 


(78a) 

(78b) 

(79) 

(80) 

W ailr 

'F3cp= 

(81) 

)] 

ailr~] 


0 


(81a) 
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Flow on Radial Plane 

For "blie special case of flow on a radial plane^ 

= n^^ = fr = fu = ° 


and eauations (77) to (8l) reduce to 


dz 


i- -VpM, 


dilr _ 




= rt’pW^ 


li =. I + i 


r/i dtf . 

“Lys^y wJ 

/ Wr2\a2,lr . (l + + 0 

\ a2/r2d(?? ^ ^ 


(82) 

(83) 

(84a) 

(84b) 

(85) 


(86) 


where 


M == 


or 




d In b’ ds* 1 


(- 


N = - - 


1 d In b' 1 Ss* 1 

^ tW ^‘2- 


dl a2-w2 ^ ds^ a2-w2 2 cd 




a r 


Wr^ ^ 




a2 ^^r 


and 


d2,|r 1 dUr 1 d2t /d In b'p dt 1_ d la b ^ , _ 


(t'p)' 
1 dT[r 
r dv 



0 


(86a) 
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Alternative Form of Equations for Flow along 
Surface of Revolution 


The equations given in the preceding sections are obtained for 
turhomachines to avoid an infinite value of the partial derivative with 
respect to z. Difficulty still exists in using either of the systems 
in the case of a mixed-flow type machine with an axial inflow and a 
radial outflow. For solutions of general surfaces, this difficulty 

can he avoided by dividing the machine at the middle of the flow path 
and using the first system at the inlet portion and the second system 
at the exit portion. If the surface can be approximated by a sur- 

face of revolution, it is convenient to use a set of orthogonal coordi- 
nates I and cp, where I is the arc length of the generating line of 
the surface of revolution in the meridional plane and cp is the usual 
cylindrical angle (fig. 8) . Because 


and 


^r dr . 

^ a — = sin a 
¥2 dZ 


dz 

^ a ^ = cos o 
¥2 dZ 


then, for use with the first system. 


9q Sq 


and, for use with the second system. 


3q „ ^q 

~ =r CSC a 
9r dZ 


(87) 


( 88 ) 


(89) 


(90) 


By use of the preceding relations, the equation of motion in the circum- 
ferential direction as given by either equation (56) or (79) for the 
two systems, respectively, becomes in both cases 


1 

r ^9 dz 



sin o 


¥, 


dl 

d9 




(79a) 


which agrees with the results obtained in references 29 and 30 in a 
different manner. The subsequent equations given in these two refer- 
ences can be modified and used for such surfaces. (The last teim on the 
left side of equation (79a) represents the rotationality of the absolute 
flow and is not included in reference 29, which is derived for Irrota- 
tional absolute flow.) 
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By compaxlng the integrating factor h used herein and the thick- 
ness of the stream filament of revolution t used in these two refer- 
ences, it is seen that the two play exactly the same role in the con- 
tinuity relation. Although h is obtained matheimtically as an inte- 
grating factor, physically it may then he visualized as the thickness 
of the stream filament in the r or z-direction for the two systems, 
respectively. The use of h herein is, of course, more general in that 
it varies two -dimensionally over the surface in the general case, where, 
as in references 29 and 30, t is considered a function of Z only. 


EQUATIONS GOVERNING FlUID FL0V7 ON Sg SURFACES 


In the preceding section, it was shown that the determination of 
the flow on Sp surfaces req^uires a knowledge of the radial var-iation 
of the velocity con^ionents . This knovrledge can be obtained by following 
the fluid motion along relative stream surfaces of the second kind, Sg- 
On S 2 J the relations (28) to (31) also hold. These relations, how- 
ever, will now be used to eliminate the Independent variable Cp; that 
is, any quantity q on S 2 is now considered as 


Accordingly, on S 2 


q = q[^r,z,(p(r,z)]] 


^ _ dq °r 1 5q 
ar ” ~ n^ r ^ 


3q _ 5q ^ 1 Sq 
3z ~ ^ ~ n^ r 5^ 


and along a streamline on S 2 


Uq 

Dt 


= Wr 


3r 




3Z 


(91) 


(92) 


Equations of Continuity and Motion 


Equations (30) and (9l) are used to change the continuity equa- 
tion ( 35 ) to 


1 3(p¥rr) a(pW2) 
r 8r 3z 


P C(r,z) 


(93) 



6V 


KAjCA TN 2604 


39 


where 






C(r,z) - - °"r . 3Wu 


^Cp 


three perpendicular directions are 





(94) 

of motion (l4a) 

in the 


MufSCVu: 

r L Sr 


SWj, 


+ W 


/aWr _ ^z\ _ 

z \5z dr / 


ai 


as 


^ ^ SF 




- w. 


^ui’) 

at^r 

TT ^Z 

awu\ 

dr 

■ acp_ 

\r d(p 

■ 

/aWr 

^z 


awA 

\az 

“ ST 

) *1 \r acp 

■ 3T/ 


_ 1 ai T as 

“ r 5^ r ^ 


> ( 95 ) 


ai 


+ T 


as 




the motion on 82 ^ equations (95) are reduced to the fol- 
lowing form by using equations ( 9 ), (I 6 ), (3l), and (91): 


Wu 9 (Vut) 


+ W, 


\8z 


3r ^ V 8z 

Wr 9 (VuT) a(V^,r) 


_ 91 n, 93 

ar / " ar ^ ar ■*■ 


ar 


-w 


r 

aw 

ar 


= F, 


az “U- 
z\ ^u ®(^u^^ 


or FyT ss 


D(V^r)' 


Dt 


az 


^ + T^ + F 
az + a^ + 


(96a) 

(96b) 

(96c) 


where F is a vector havxng the unit of force per unit mass of gas 
defined by: 


F = - (1^ . T 


1 1 ap 

ny^r p ^ " 


(97) 


A similar result is obtained for the equation of motion ^ n the form of 

( 2 ) : 


W 


aw^ 


aw^ 


+ w„ 


r ar "z az 

^r 9(V^r) Wg a(Vyr) 


lap , ^ 

r " p ar ■*■ ^r 


ar 


aWz 

w„ :^ + W- — 

^ az ar 


r az 
aWz 


■•U 


= - - ^ + F 
p az z 


J 


( 98 ) 



40 


mCA W 2604 


Because the vector 


By the use of eq^uations 
on an Sg surface that 


F is normal to the 

^r^r ^u^u ~ 

(99) and (96), it can 

Dt Dt 


S 2 surface, 

0 (99) 

he shOTm for steady flow 

(99a) 


This result is the same as that obtained for the surface. Again, 

for the present problem of steady relative flo^r on the 83 surface, the 
relation (99a) can be taken either as one of the equations of motion or 
to represent the relation (99) . In other words, there are only four 
Independent relations among equations (96a), (96b), (96c), (99), and 
(99a) . In the foUoiri-ng development, it is found convenient to use 
equations of motion in the form of equation (96) , not only because 
hx/br is zero in many design problems (whereas ^ O) , but also 

because equation (96) leads to a form capable of a rigorous solution 
for both subsonic and stq)ersonic flow and shows clearly how the various 
design factors affect the three-dimensional motion in general. (See 
equations (IO 6 ) to (114) that follow.) 


In a manner analogous to the surface, the continuity equa- 

tion ( 93 ) is put into the form 


a(rBpWj,) 

ar 


+ 


a(rBpW^) 

az 


= 0 


( 100 ) 


by the use of an integrating factor B, which is related to C by the 
foUowing equation: 


or 


DlnB ainB.^^aiJiB 

Dt ^ ar z; 32 


( 101 ) 



(101a) 


Equation (lOO) is the necessary and sufficient condition that a 
stream function ilr exist and 


ar 


= rBpW^ 


(I02a) 
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az 


- rBpWj, 


(102b) 


The difference in ¥ at two points j and k on the Sg s\arface is 




rBp(Wz dr - Wr dz) 


Similar to' the flow on the S-j_ surface, the preceding equation 
indicates that B is proportional to the angular thickness of a thin 
stream sheet whose mean surface is the stream surface Sg considered 
herein and whose variable circumferential thickness is equal to rB. 
Indeed, if the mass flow into and out of the element of such a stream 
sheet (cut between two planes normal to the z-axis, and a distance dz 
apart and between two cylindrical surfaces dr apart (fig. 7(c))) is 
equated to zero and the distances dr and dz approach zero as a limit, 
the following equation is obtained; 


3(TpWr) a(TpWz) ^ 
3r ^ 3z 


(100a) 


Con^jaring this equation with equation (lOO) and considering the mass 
flow relations show T to be proportional to rB. This proportionality 
means that B can be physically intepreted as a quantity which is pro- 
portional to the singular thickness of a stream sheet whose mean surface 
is the S 2 surface considered herein. With this interpretation, B is 
immediately seen to be closely related to the angular distance between 
two neighboring blades. In actual calculation, only the ratio rB to 
(rB)^ or T to Ti is in 5 )ortant, and it is also easier to obtain the 
variation in rB from the distance between adjacent streaml in es 
obtained on S]_ surfaces than to evaluate B/B;j^ by equations (lOla) 
and (94) using data obtained on Sp surfaces. 


Principal Equation for Case with V^r Given 

In the solution of flow on an Sg surface, the continuity equations 
and the equation of motion in the radial direction are combined to form 
the principal equation. The principal equation will now be obtained for 
two DKiin groups of present designs in which a certain desirable varia- 
tion of the angular momentum of the fluid V^r and of the ratio of 
relative tangential and axial velocity are prescribed on the 

face, respectively. These equations can also be used for the solution 
of a direct problem, in which the same information obtained on Sp 
solutions of a previous cycle is used as known values in the Sg 
solution. 
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For the first group, the foUo\ri.ng eq^uation is considered knorira: 

VyT = G(r,z) (103) 

Among this group of designs axe the free -vortex design (in which G is 
simply a function of z), the more general "solid-body rotation" design, 
the "symmetrical velocity diagram at a3J. radii" design, and others (foi 
exaii 5 )le, see references 17 and 18) . 

From equations (102) and (45), 


3^2 a^Tjr _ / 1 1 ah 

Bp r " ^2 ar 


3s* 3 In B\ 3¥ 


ar ^2 
3r 


3r 


3r 


\ 

j ar 


3% , / 1 3h 3s* 3 In B\ 3f 

^ az ^ 2 2 az ^ 3z ~ az ) az 

az \ a / 


But from equations (9) and (lOl), 


1 / \ —2 

H - I + ^ - I (rBp) 


3r j 


^ mf 


(104a) 

(104b) 


(105) 


Differentiating with respect to r and z gives 

S = i (- • <».W> (^ - r ) - 

(rBp) ^ 

Substituting the preceding equations into equations (104) and adding 
give 


[a^- .BP - g) = (a^ V) 


a^Tir a^ nf 

y + 

ar 


3r 


^+r 

a2 

a W^2^2r2'^ 

- a2 

/'a lx B 3 s*M 


r 

ar V 2 / 


\ 3r 

3r /J 

r 3 



In B 

as*\ 1 

3HT 

L 3z 


- 2 

3z 

-wj\ 

3z 


( 106 ) 
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Substituting equation (l06) into equation (96a) and dividing "by 
a^ yield the following principal equation for the fluid flo\7 on sur- 
face Sg: 


('-¥) 9 -= 


^r^z 3^'f' 


qZ ax^z 


(>■¥) 


3 z 2 3r 3z 


where 


3 In B , 3s 

M h — 

3z 3z 


N= - i - + 

r 3r ar 


a2-(Wr2+W22) 


2 


-i/ii _w 

a2 \9z 3z J 
a2 \3r 3r 


+ ca‘^r + 


(107) 


From the coefficients of the second derivatives^ the principal equation 
is seen to be hyperbolic or elliptic when the meridional velocity 

greater or less than the speed of sound, respec- 
tively. For the elliptic case, it is again convenient to write the 
principal equation in a slightly different form. From equation ( 101) , 


rBp — = — - 1 ^ _ a In Bp S¥ 
^ 3r gj,2 r 3r 3r 3r 

aWj, 32 ^ a In Bp 

^ 3z ^2 az az 

az 


J 


Substituting into equation (96a) results in 

a2y 1 / a In Bp ^ ^ a In Bp 

3 j ,2 “ r ar ^ ar 3r 3z az J 


(rBp) 2 p^u 

3J 

ar 


ar 


+ T 


ar 


= 0 


(108) 


(I07a) 
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Witli the vaxiation of or prescribed by the designer in an 

inverse prohlem or taken from the previous calculation in a direct 

problem^ the meridional velocity con^onents are determined by eq.ua- 
tions (107) and (l07a) . (Other equations are used to detennine various 
terms involved in the coefficients M and N.) 


Principal Equation for Case -with Given 


In the second group of designs, the follojdjig relation is pre- 
scribed on an S 2 ^in surface (for exanple, see references 17 and 18) ; 


Wu 

Wz 


g(r^z) 


(109) 


In order to result in blades with the mean blade surface composed of all 
radial elements (for high-speed rotation) , it may be desirable to spec- 
ify a mean Sg surface consisting of all radial elements. Then 


¥ 


u 



r 


Sl(z) 


(110) 


Similarly, in order to obtain a cooled turbine rotor blade with minimum 
twist, the follOTTing function may be specified on S 2 ^ui: 


Wu 

Wz 


g2(z) 


( 111 ) 


In application to direct problems, one of the preceding relations is 
obtained from the S]_ solution in the previous cycle and is considered 
as given in the S 2 solution. In both inverse and direct problems, 
I'Tlth the relation between ¥^ and ¥g given by these equations, all 

three velocity con^onents axe to be combined into the main terms of the 
principal equation as follcnre: Substituting relation (109) into equa- 

tion (96a) gives 


(1+g^) 


3¥„ 


ar 



+ 2cDg + 


1 


+ I + F j 
ar ar y 


( 112 ) 


InsteadL of equation (lOS) , 

- I (rBp)-=[(l+g2) 


I + 



2 


(113) 
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shoxild now "be written. Differentiating with respect to r and z, 
ccmibining with equation (104), and substituting into equation (U2) 
give the following form of the principal equation: 


(l+g2) 



^ -2(1+82) 
dr^ 


WrWz 

a2 3raz ^2 J 


?^h.n|I+mP=o 


3z' 


ar 


3z 


(114) 


where 


3 iQ B ^ a^\ 

3z 3z *" 0^2 \3z 0^2 ^ 3z J 



-L 

a2 Ur 


ar 


(O^r 


3s 

ar 


^ + F 



+ 


This equation is hyperbolic when the relative velocity is supersonic, 
elliptic when the relative velocity is subsonic. For the subsonic case, 
a form of this equation more convenient for ccarrputation is obtained by 
substituting equation (108) into (ll2): 



It may be noted that for both groups , equation ( 96a) ■ rather 
than (96c) is chosen to obtain the principal equation of the present 
problem, because F^. is always much smaller than F^ in axial 
machines and F^ is zero or nearly zero on Sg surfaces for high-speed 
centrifugal and mixed-flow impellers whose mean blade surfaces are 
usually composed of all radial elements. (For low-speed centrifugal 
impellers, equation (96c) can be used to form the principal equation 
in a similar manner.) 
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Procedure of Solution 

Althou^ the equation of motion (96a) is chosen to form the prin- 
cipal equation, other equations are to he used to obtain the various 
terms involved in the principal equation. As in the case of general 
Sq s\irfaces, there are ten basic variables to define the flow and the 
shape of the S 2 surface. They are: B, F^., F^, F^^ s, 

and I (or p). B is considered given, (in the direct problem, B 
is evaluated directly from the distances between adjacent streamlines 
or according to equation (lOOa) using the value of C obtained on Sq 
surfaces; in the inverse or design problem, B is estimated (refer- 
ences 29 and 35) from the blade thiclmess as desired from blade stress 
and other considerations . ) On the other hand, there are seven inde- 
pendent relations in one energy eqioation (21); three equations of 
motion, one of equations (107), (l07a), (U4), or (U4a), and equations 
(96b) and (96c); the orthogonality relation between T7 and F, equa- 
tion (99a); and the two equations relating f and velocity, (l02a) 
and (l02b). 

Direct problem . - In the direct problem, two alternative procedures 
may be used. If the shape of the Sg surface (detennlned from the 
data obtained on Sq siarfaces) is considered as given in the present 
S 2 solution, two additional relations between the n- or F-components 
completely define the problem. The procedure of calculation is as 
follows; 

( 1 ) Use equations (20) and (21) to determine the variation of 
s and I. 


( 2 ) Compute 


from the orthogonality relation as follows: 


= - ( — + — U 


( 3 ) Compute Fy^ from equation (96b). 

( 4 ) Solve the principal equation. 

( 5 ) Obtain ¥j. and ¥g from equations (l02). 


If only the tangential velocity or the relation (109) is taken 
from the Sq solutions of the previous cycle and is considered as 
given in the present S 2 solution, one more relation is available 
between the F-components such as that which exists between the 
f-components on the Sq surface: 


F • V X F = 0 


(115) 
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Writing eqimtion (U5) in scalar form and using the relations (3l) 
and (9l) give 


9 ^ pr \ 

9r \FuV 8z \Fui*/ 


(U5a) 


By integrating along a constant z-line, equation (U5a) provides the 
following relation to determine the value of to he used in the 

principal equation from the values of Fu and Fz: 



If Fj. = 0^ at Zq 



The procedure of calculation is as follows: 


(U5b) 


(H5c) 


(l) Use equations (21) and (99a) to determine the variation of s 
and I. 


(2) Compute Fu and F 2 from equations (96b) and (96c). 

( 3 ) Compute Fj, from equation (U5b) or (U5c). 

( 4 ) Solve ¥ from the principal equation. 

( 5 ) Compute Wp and Wg from equations (l02a) and (l02b). 

Inverse problem . - In the Inverse or design problem of a finite 
number of thick blades^ in addition to the blade -thickness distribution 
or its equivalent either equation (l03) or (l09) is prescribed on 
a mean stream surface S2^jn* appear that still another rela- 

tion can be prescribed on the mean surface. The differentials of the 
coordinates of the siarface are now governed by 


Fj, dr + F^r dqj + Fg dz = 0 


(116) 
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however, and, in order that this differential equation will lead to an 
integral surface of the form repres^ted hy equation (28), F must 
satisfy the condition of integrahility as given hy equation (115) 
(reference 34). An ejcpression similar to equation (U5a) for the case 
of an infinite number of blades was first pointed out hy Bauersfeld 
(reference 2) in a discussion of the Lorenz paper (reference l) . In 
effect, it restricts the freedom that the designer has in prescrihing 
the velocity components of the flid^d on the surface. Hence, in the 
Inverse problem of a finite nunlier of thick blades, in addition to the 
blade thickness distribution or its equivalent B, the designer can 
specify only one relation on the mean stream surface, which relation may 
be either the tangential velocity as given by equation (103), the flow 
angle between the tangential and a-xi.al velocity as given by equa- 
tion (l09), the axial velocity, or any other reasonable relation that 
will lead to a solution of the set of equations. 

In the preceding consideration, the hub and casing shapes are also 
prescribed by the designer in the inverse problem. Alternatively, the 
prescription of the hub shape can be replaced by a prescription of 
another relation at the casing, thereby fixing the shape of and the flow 
along S2,m "tlie casing entirely. The flow is then extended to the 
hub and the last streamline gives the hub contotur (reference 19). 


i^roximatlons Involved in Through-Flow Theory 

“When the equations previously derived in reference 18 for a large 
number of thin blades are compared with the corresponding equations 
derived herein along a stream surface, the two are obviously exactly 
the same if the ordinary derivatives used in reference 18 are replaced 
by the present partial derivatives following the stream surface, and 
if B is equal to 1 or if the VEirlation of B along the flow path 
is zero. In the interpretation of the throu^-flow solutions as the 
flow SLlong a mean stream surface (which divides mass flow into two eqxaal 
parts clrcumferentieQly) or as the flow along the' mean channel surface 
(geometrical mean), the first difference can easily be removed by simply 
interpreting the values obtained in the solution as those along the stir- 
face rather than in the meridional plane. The second condition, however, 
is satisfied only when the circumferential variation of all the velocity 
ccmponents approaches zero, or when the circumferential derivative of 
the tangential velocity and the ratios of n^ and nz to n^ approach 
zero (see equation (94)). 

Besides the use as a limiting solution in general and to give cer- 
tain trends where the contribution due to the finite number of blaides 
is sman. or constant, the throu^-flow calculation should be properly 
modified by the factor B in its application to actual turbomachines 
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of a finite number of thick "blades. As B can "be physically inter- 
preted as the ratio of the local angular thickness of the stream sheet 
to its inlet value^ a good approximate value can he obtained hy solving 
the two-dimensional flows on a number of. stream surfaces of revolution 
starting at different radii at the inlet. For. the subsonic flow in the 
turbine cascade reported in reference 29 and for the supersondc flow in 
t\70 impulse bladings investigated in reference 30^ the reductions in 
angular thickness from the inlet value along the mean streamline are 
seen to be a chordwise average of 4 and 9 percent more than the reduc- 
tion in the channel width, respectively. Also, in the subsonic case, 
the Influence is extended a certain distance outside the blade row. The 
incluaion of this factor B, even if it is approximate, should give a 
much closer answer than that obtained with B taken as 1. 

In this interpretation of the infinite number of blades solution 
as the solution of through flow along a particular stream surface 
between two adjacent blades, the distributed "body force" F has a 
definite meaning, as given by equation (97). (For an infinite number 
of blades, F becomes the blade force.) For blades with large turning 
and large radial twist, as in a free -vortex turbine, the Influence of 
the radial consonant of F on the flow is not negligible. 

CIRCDMFEREMTIAL VAEIATION OF FLUID PROEEKTIES BI USE OF POWER SERIES 

In general, the blade -to -blade variations of fl\ild properties are 
to be obtained from calculations on Sq surfaces. "When the twist of 
the Sq siarface is large, some other method of obtaining the blade-to- 
blade information is desirable, For subsonic irrotational absolute 
flow, this Information can be obtained by extending the solution 
obtained on the mean stream surface in the circumferential direction by 
the use of power series (without the consideration of the shape of the 
Sq flow surfaces) . The various derivatives involved in the series are 
obtained from the flow condition on the mean stream surface. The hi^er 
the solidity and the thinner the bleide sections, the fewer are the terms 
required for a given accuracy. Results obtained in references 29 and 36 
indicate that only three terms in the series will be required to give 
sufficient accuracy for high-solidity turbines and centrifugal, com- 
pressors. 

The series method will also be \zsed in one of the two methods of 
the inverse solution in which the flow obtained on the mean stream sur- 
face is extended out circumferentially. 

Denoting the absolute vorticity VXV by | and using the rela- 
tions (16), (9l), and (97) in equations (l9) give 



50 


MCA THf 2604 


er 

1 

1 

r^z ^('^'ur)"^ 

X 1 a(Vur) 

Fz 

^u 

(117a) 



dz 

~ r 3cp “ r 

Fur 


eu 


aHj. 

aWz Fz 5Wp Fp 

SWz 


(117b) 

= " 5?"; 

3z 

" ar **" F^^r Sep " F^r 

Sq5 


5z 

1 

r S(V^r) 


1 a(V^r) Fr 

1 

SWr 

(117c) 

~ r 

_ ^r 

-w] 

~ r ar F^^^r Scp 

r 

Sep 

From the preceding equations. 






3^2 

Frir + 

S'ulu 

+ Fziz Fz a(V^r) 

F 

■^r 

a(v^r) 

(118) 

az 

ar 


Fu 

" F^r ar 

Fur 

82 


This equation means that the apparent vorticity, which is ohtainei hy 
differentiating the velocity on the mean stream surface with respect to 
the coordinates, is not zero even if the absolute vorticity is zero or 
tangent to the mean stream surface. Substituting equations (U7a) 
and (U7c) into equation (94) results in 


r 5 ^ ■ P~L^u- 


a(V^r) a(T^r) . _ F^ ^ F^ ^ 

r 3r F^r 3z ^ * Fu " F^ ^ 


( 119 ) 


Substituting equation (ll9) into equations (U7a) and (I17c) gives 
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The first derivatives of h or p, and p can be obtained as follows: 
From equation (97), 

r 5^ r 5^ pr ^ 


( 122 ) 



MCA TW 2604 


51 


Or from equation. (9), 
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can he obtained by using equation (l 2 b): 


(I22a) 
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The second derivatives of the fltiid properties with respect to qp 
can be obtained in a similar manner. Differentiating the continuity 
equation (l) with respect to cp and dividing by r give 


J .2 ~5r~5^ ^ J .2 ^ r dz dcp 

Equations (9l) are used to change equation (l24) to 


(124) 
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Differentiating equations (U7a) and (ll7c) with respect to cp and 
dividing by r result in 


i ^ ^ 1 aSr 1 

p 2 r ^ ^2 ^z 3cp 


1 ^^r _ l^^z 1 ^^(V^r) 
r 2 ^ " r ^ ^2 3r dcp 


1 i .a. 1 
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Substituting equations (l25) and (126) into (l24a) and noting that F 
is perpendicular to -W give 

~^r J_ _a ^(yu^) + 

j,^^2 3r §9 F^r az pj.2 ar 

1 a ^(P^z) 2C a to p 

pr az 5^ r ^ F^r ^cp Fy^r ^ 

Equation (l03) is to he used in equaticms (lOl) and (l02) to obtain the 
second derivatives of Mx and The second derivatives of h and 

p are obtained from equations (9) and (123) as 


(127) 



2l. ^ i ^ 

r2 acp? x^ Scf? 2r^ ^ 


(128) 


i P ^ 1 h _ ^ 

Scf? a^r^ Bcp? r^ 3cp^ 


(129) 


Similar formulas can be obtatoed for higher-order derivatives. At 
a fixed value of r and z, the velocity components, h, and p at a 
short angular distance away from the mean stream surface Sg can then 
be obtatoed by a Taylor series: 

(cp-cp )2 (Cp-Cp 

l(cp) = q(opJ + (^-%) ^ ^ <l'"('Pm) + • • • 

(130) 

An alternative way to obtain density is to use equation (l45) (to be 
given subsequently) after the other fltiid properties are determined. 
Obviously, the preceding equations are most useful when the flow is 
isentropic with vorticity equal to zero. Otherwise, the variation of 
vorticity along the mean stream surface has to be determined first. 

At present no such method is available. It appears, however, that the 
method of Squire and ¥toter (reference 37) and Hawthorne (reference 38) 
may be generalized to compressible flow for the variation of vorticity 
along a mean stream surface to torbomachtoes. 
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STEPS FOR CCMPIETE SaLOTICOSB OP THREE-DIMEHSiaHAL 
DIRECT AHD DIVERSE PROBEaiS 

Di general, the solution of the three-dimensional direct and inverse 
problem involves the nee of both Sp and S 2 surfaces. In the direct 
problem, starting with assumed flow surface, the solution is obtained 
throu^ the successive (alternate) use of the two kinds of flow sur- 
face, although a satisfactory approximate solution may be obtained in 
one or two complete cycles. The use of an approximate method of solu- 
tion to get a good starting value on each surface will shorten the length 
of computation. For inverse problems, the process is usually shorter. 

The calculation will start on the S 2 m surface on which either a con- 
dition on the fluid velocity or the shape itself is prescribed and an 
estimated value of B for a desirable blade thickness distribution is 
used. After the solution on the S 2 ^m surface and its shape are obtained, 
the blade coordinates are obtained by extending the solution circumferen- 
tially either by the series method or by the method given in reference 35 
using the variation of the distances between the streamline obtained in the 
S2,m surface. Because it is in 5 )ortant only to obtain the right order of 
magnitude and the rl^t kind of variation (three-dimensionally) of the 
blade thickness, the first solution may give satisfactory results. The 
velocity distribution on the blade surface is controlled directly by the 
one relation specified on the S 2 jm surface and the variation of B. 

Suitable procedure is subsequently suggested for the solutions of 
direct and inverse problems with either irrotational or rotational inlet 
absolute motion, at design or off -design flow conditions, for tiurbo- 
machines having various wall configurations (fig. 6) . 


Direct Problem 

Axial turbomachines with nontapered straight walls . - In this type 
of machine, it is desirable to start the ccmputation on Sp surfaces, 
because with short axial blade length, the total deviation of the 
surface from the cylindrlcELl surface is relatively small, especially 
along the hub and casing walls. 

The following steps are therefore suggested: 

(a) In the initial calculation, the flow siirfaces are assumed to be 
cylindrical and the set of equations (60) to (65) derived for cylindrical 
flow or the approximate method given in reference 39 can also be used to 
obtain the streamlines and circumferential variation of fluid state on 
Si surfaces at three or more radii. 
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(Td) From the data obtained in step (a), an Sg stream surface about 
midway between two blades is constructed by connecting the streamlines 
which divide mass flow on the Sp surfaces in the same percentages. The 
direction n-umbers of the surface and the ¥u. and Wg at the surface 
obtained in step (a) give the starting value of by use of equa- 

tion (31). The factor B is evaluated either directly from the angular 
distances between streamlines obtained in (a) or according to eqm- 
tion (lOOa) with C evaluated from the information obtained in (a). Its 
value at other radii is obtained by interpolation or by proportioning 
according to the channel-width ratio. Calculation of the flow on this 
surface is then made by the use of equations (9l) to (U5). For subsonic 
flow with Irrotatlonal inlet flow, the solution obtained on the 
surface is easily extended circumferentially by series expansion using 
equations (U7) to (130). The values obtained can be further adjusted 
to fit the given blade (reference 39) and can be used in a more accurate 
second calculation on Sp surfaces in the next step. For subsonic flow 
with large rotationality at the inlet and supersonic flow with signifi- 
cant check caused by the blade entrance angle, it is more desirable to 
obtain the information on circumferential variations by the use of two 
or more Sg surfaces at or near the two blade surfaces. 

(c) The radial variation of fluid state camputed from the solution 
obtained in step (b) or the variation of the radial distance between 
streamlines is used to determine the factor b and used in the principal 
equation (48) for a more accurate determination of Sp surfaces and the 
flows thereon. The general equations (32) to (5l) should now be used 
for the Sp s\arfaces located between hub and casing, if not at or near 
these walls. 

(d) The calculation of Sg surfaces can again be repeated and so 
forth. 

If the inlet flow is quite rotational, so that the Sp surfaces 
along the walls and the Sg stirface near the blades may turn around the 
comers, these surfaces should be chosen at a short distance from these 
boundary trails as shown in figure 3. By the use of these two vindR of 
surface, the secondary flow caused by a rotational inlet profile or by 
the tu rning of the blades is included in the complete solution. 

Axial turbomachines with tapered or ctirved walls . - The steps 
involved here eire quite similar to those of the preceding case, except 
that for the i ni tial calculation of Sp surfaces along or near the 
tapered or curved wall, either equations (52) to (59) are to be used, or 
equations (13* ) to (23’ ) given in reference (35) can be used in the 
manner given in reference ( 39 ) . 
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Badlal- and mlxed-flovr type -tnirbcanachlnes with curved -wa-lls . - In 
this type of machine it is not desirable to start the computation on the 
Si surfaces because the flaw surfaces near the walls may deviate con- 
siderably from surfaces of revolution because of the long flow path. On 
the other hand, the solidity of the blade is very hl^ and the blade 
section is uniformly thin. As a resiilt, the shapes of the S2 surfaces 
are closely related to the blade shape and the factor B can be esti- 
mated relatively accurately from the blade thickness distribution. The 
following steps are therefore suggested: 

(a) The computation is begun on the S2^m surface. For subsonic 
irrotatlonal inlet flow, computation need be made only on a mean 83 
surface and the solution can be extended out circumferentially by equa- 
tions (117) to (130). The approximate method given in reference 40 can 
also be used in the initial calculation. For subsonic flow with rota- 
tional inlet profile and for supersonic flow it may again be more 
desirable to compute two or more 83 surfaces between the blades. 

(b) The data obtained in step (a) may be used to make calculations 
for three or more Sp surfaces between hub and casing walls. 

(c) The processes (a) and (b) can be repeated until the desired 
accuracy is reached. 


Inverse Problem 

Conditions prescribed on mean stream surface . - In the inverse or 
design problem it is most convenient to consider a mean stream surface 
of the 82 kind about midway between two neighboring blades to be 
designed (figs. 3 to 5). From the results developed previously for 
such surfaces, it is seen that in addition to the factor B, the 
designer can specify only one relation among the fluid properties on 
that surface, which can be either a velocity component, a relation 
between two velocity components, or one other reasonable condition. 

The factor B essentially controls the blade thickness distribution, 
whereas the relation specified on the surface essentially controls the 
mean camber surface of the blade. From a consideration of strength and 
Mach number in general, and the requirement of coolant passage in the 
case of cooled turbine blades, the designer always has a very good idea 
of what kind of blade-thickness distribution he wants. ¥ith this thick- 
ness distribution, the ratio of pitch minus circumferential thickness 
of blade to pitch can be obtained. After correcting this ratio with 
seme known relations between this ratio and B (such as those given in 
references 29, 30, and 35), especially near the leading and trailing 
edges, it can be taken as the factor B in equation (lOl) . Then from 
the type of velocity diagram or a certain feature of blade shape 
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desired, a relaticai along tlie mean stream surface S 2 jji can te pre- 
scribed and coord in ates of the mean surface and the flow on that surface 
can he solved at the same time hy equations (lOl) to (U5) . It may he 
noted that in this process, the designer still has, in general, a little 
freedom in choosing the value of Zq in. equation (U5c) . For a rota- 
ting hlade, zq is usually taken somewhere near the center of gravity 
of the hlade section, whereas for the stationary hlade, the position of 
Zq can he utilized to control the magnitude and dlstrihution of Fj. 
in the most favorable manner. 

Boundary conditions for mean stream surface . - In the solution of 
this S 2 ,m surface, the boundary conditions are a little different for 
subsonic and supersonic flow. For subsonic flow, not only the varia- 
tions of the stream function at stations far upstream and downstream 
are given, the meridional contours of the huh and casing walls are also 
given (these contours can he determined hy approximate calcialations 
from hlade row to hla.de row such as given in references 17 and 41) . For 
supersonic flow, the variation of the stream function and its nomnal 
derivative is prescribed on an initial curve, which is not a character- 
istic cu 2 Tve. Then either the huh and casing contours are prescribed, 
or only the casing contour hut with one more velocity component along 
the casing is prescribed. In the second case, the flow is extended 
toward the axis of the machine and the huh contour is determined hy the 
shape of the last streamline for the required mass flow. 

Determination of hlade shape . - For siibsonic irrotational flow, 
the solution obtained on the mean stream surface can he extended out 
circ\ 3 mferentially hy using equations (117) to (130). The hlade sur- 
face can he then determined as foUo^fs: 

(a) The position of the mean stre^ surface is first determined 
hy solving the circumferential coordinate as a function of the axial 
coordinate at several radii. With the circumferential coordinate 
measured from the radial element of the svurface chosen at Zq, equa- 
tion (II 6 ) gives, at a constant r; 


mn. - (^)z=Zo 


Jzo 



(131) 


(h) The hlade coordinates (r, cp ) will first he chosen at one 
station Zq as follows (see fig. 9): The mass flow passing through 

the Zq plane between the mean stream surface and the tentative suction 
or pressure surface is coD 5 )uted as follows: 
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r^c 

^ = L p W^r dr dcp (132) 


Mp 





p Wz r dr dtp 


(133) 


Because of the Inaccuracy in B for the blade -thickness effect, the 
mass flow obtained will be a little different from that required. The 
blade coordinates cpg and cpp as functions of Zq and r are modi- 
fied until the mass flow checks. It is not important that Mg and 
are a little different from one-half the required mass flow as long as 
their sum is equal to the total mass flow, but once the division is 
chosen, it should be maintained at other z- stations. 


(c) The blade coor din ates obtained at z = Zq are extended 
upstream £ind downstream according to the velocity components evaluated 
at the surface. For example, for a short distance z - Zq away, the 
changes in the blade surface coordinates r and cp are 


r = ro + (z - zo) ' (134) 


- '’’o (f - -o) (1=5) 

After r and cp are thus obtained, the total mass flow may be checked 
again by equations (l32) and (l33). 

When the blade coordinates are obtained close to the leading and 
trailing edges, they can be faired in according to some standard shapes. 
A blade shape is therefore obtained in which the three-dimensional flow 
of the fluid is considered. The right kind of three-dimensional blade- 
thickness distribution is obtained and a good knowledge about the flow 
on the blade surface is also ayailable at the same time. The’ data 
obtained in the solution can also be used directly for a more accui*ate 
and detailed determination of the velocity variation around the nose of 
the blade, for a relatively quick check of the series approximation, or 
for improvement, if necessary, of the inverse solution throu^out by 
the method given earlier for solving the direct problem. This process 
seems to be the quickest way of establishing some standard three- 
dimensional flow variations for typical designs of blades from which a 
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good approximate method for routine design calculations can he estab- 
lished and of providing a basis on which the viscous flow along the 
blade surfaces and hub and casing walls can be analyzed. The results 
given in references 29 and 35 Indicate that for blades of high 
solidity, three terms in the series give sufficiently accurate 
results . 

For subsonic flow with vorticity, the clrc\imferentlal extension 
cannot be accurately made at the present because of lack of adequate 
methods for the determination of vorticity variation along the mean 
stream surface (Sg^m) • ^ estimate of this variation can be made, 

however, and the solution can be cheeked later. An alternative method 
is to use the shapes of the streamlines and the distances between them 
obtained in the S2^m solution and to design the blades with the 
assun 5 )tion that the flow surfaces are surfaces of revolution by the 
method given in reference 35. Inasmuch as the rotational! ty of inlet 
flow is usually serious only in later stages of a multistage compressor 
where the hub-tip radius ratio is high^ this assumption is reasonably 
good. 


For supersonic flow, the flow in the mean stream surface S2^m is 
also determined first. If the shock due to the entrance wedge ari^e 
is small, an approximate solution of the blade shape can also be 
obtained by .the series method neglecting the finite jump across the 
shock or using an estimated value. The improvement of the flow varia- 
tion for the resultant blade is then more important than that in the 
subsonic case. Local modification of the blade shape can also be made 
if the velocity distribution on the blade obtained is unsatisfactory. 

A better method is to use the shape of the streamlines and the dis- 
tances between them obtained in the S2,m solution and to design the 
blades assuming flow surfaces of revolution according to the method 
given in reference 30. 

The processes described here for the three-dimensional solution 
have been and are now being used to analyze the compressible flow 
through a numiber of compressors and biirbines. Seme of the results 
obtained are given in reference 42. 


GENERAL METHODS OF SOLVING ERINCIPAL EQUATION 

In the solution of the Sp surface for the direct problem and of 
the S 2 surface for both the direct and Inverse problems, the mairi 
calculation is the solution of the principal equation, which is a 
second-oixLer, nonlinear paxtial differential equation in t\TO independent 
variables. The case when the principal equation is elliptic will be 
considered first. 
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KLULptic Case 

A common fonn of the principal equation is -written as follows: 


nan a^2 a^ 


a^Tif 
^ T 4_ 
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an 


In hp 


+ K 


In 




MjLhp 


P 


M2(hp)2 = 0 

;pi equation (l36), >Jr and h are used for "both and Sg surfaces; 

n denotes <p for the Si surface and r for the S 2 s-urface; 

^ denotes z or r for the Si surface, and z for the S 2 sur- 
face. The "values of t]^ and K for each indi-vldual case are gi-ven 

in the following table: 


Case 

Surface 

Coordinates 

IK| 

K 

L 

Equation 



n 


ill 



1 

Si (general) 

cp 

z 

0 

r2 

“u 

Urr 

r\ 

(48a) 
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Si (surface of 
revolution) 
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2 i 

1+X2 

(59a) 
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U 

3 

Si (cylindrical 
surface) 


z 

0 


0 

(65a) 

4 

S]^ (general) 


r 

0 

1 

r2 

0 

(75a) 

5 

Si (surface of 
revolution) 

CP 


1 

1 +X^ 


(81a) 
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Si (radieil 
plane) 

CP 

r 
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0 

(86a) 
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0 
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(l07a) 
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The equation is nonlinear even in the case of incompressihle flow. In 
the numerical cor^jutation, it is convenient to rewrite the equation in 
the following form: 


where 



3\jf 


+ K 


8^ 

3Ti^ 


+ L 


= N 


a In hp ai{f a In hp aijf 

ari aq a^ a^ 


M^bp 


M2(bp)2 

3^1 


(137) 


(138) 


and is evaluated from any approximate solution at the start of the cal- 
culation and from the values of ijr and p obtained in the previous 
cycle during the calciilatlon. For simple boundary shapes for an Sg 
surface and sin 5 )le functions of J, K, and L, it is possible to find 
a Green’s function G(ti,-^, 'X., y) with its proper characteristics so 
that the solution of the problem can be written in the following form 
(for example, see reference 10) ; 


t (n^^) 




G (t);. X, y) N (x, y) dx dy 


(139) 


If the boundary wall is arbitrarily curved, it is necessary in this 
method to use the technique of conformal transformation to render the 
given boundary into a sinqiler one, such as cylindrical. Because this 
process will involve a numerical solution of the Laplace equation with 
the given boundary shape, it may be better to solve the given equa- 
tion (137) directly with the given shape by the numerical method. 
Fiarthermore, this method will be the only choice in the general case 
where J and K are complicated functions, which mahes the task of 
obtaining the proper Green’s function very difficult if not impossible. 

Finite-difference form of principal equation . - In order to solve 
the given equation (137) directly, the general numerical differentiation 
foarmula for first and second derivatives with the function value given at 
unequally spaced grid points using second- and higher-degree polynomial 
representation as given by reference 26 is used to give the finite- 
difference expressions conveniently and accurately at the grid point 
near the curved boundary. If the value of any quantity q on the 
stream surface under investigation corresponding to a number of values 
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of one of the independent variables x not equally spaced, denoted by 

xl, . . . x^, . . . x^, is given, the derivative of q (on 

the surface) "with respect to x when x = x^ may be written 


(rFq) i = y 

x=x^ Z_y n j n 

j=0 


(140) 


IThe differentiation coefficients B and the coefficients of the deriva- 
tives in the first or second remainder term have been explicitly 
expressed in reference 26 in terms of the spacings between the successive 
grid points by using polynomials of the second, third., and fourth degree 
for general nonimiform spacing throu^out and for the speclsil case neeir 
a tapered or curved boundary where only the first or last spacing is 
different from the others. For the special case, these coefficients 
have also been computed for different ratios of the distance between the 
boundaiy and the nearest point and the other spacing, from 0.1 to 1.29 
in intervals of 0.01, and are given in reference 26. For spacing ratios 
lying between these tabulated intei'vals, B can be obtained from the 
values tabulated by applying inteipolation formulas given in refer- 
ence 43, or by the direct use of the formulas. Differentiation- coeffi- 
cients B for equal intervals using various degrees of polynomials are 
given earlier by Bickley in reference 44. 


In the present fluid-flow problems, a large region must be covered 
in order to get to the boundary conditions which are always given at 
stations far upstream and downstream of the blade row. In order to 
reduce the labor of confutation, it is desirable to attempt to reduce 
the number of grid points req-ulred for a given accuracy by using a 
degree of polynomial higher than the customary second. A study of the 
expression of the remainder terms (see reference 26) and actual experi- 
ence in the present problem show that, in most cases, the use of the 
fourth-order polynomial will reduce the necessary number of grid points 
to less than one-fourth that required by the second-order polynomial. 
Wear the leading and trailing edges of the flow on surfaces of the Sq 
kind, the variation of ijf is such that accuracy is obtained most 
effectively by using small spacirg there. In such case, the grid 
pattern should be chosen at these regions first, and either be kept 
constant or be continuously increased toward the inlet and exit stations. 


With the grid pattern aind the order of polynomial representation 
selected, the coefficients B at each point can be obtained from refer- 
ences 44 and 26 for equally and unequally spaced points. Then the dif- 
ferentieil equation (137 ) at any grid point whose ijr value is T]f^ is 
replaced by the following algebraic equation: 
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+ Ki lAi) ¥ + (?Ai + Ia^) ilr^ - = 0 (l4l) 

j=0 ^ ^ ^ ^ ^ ^ ^ 

where ijr^ and ili^ denote the values of ijr on the surface considered 
corresponding to the grid points along constant ^ and constant t] 
lines, respectively. (See figs. 10 and 11.) It should he noted that, 
in accordance with the definition of the special partial derivatives, 
ijr values are those on the surface Sj whereas the grid spaclngs 
involved are just the distances along the ti- and ^-coordinates. 

Boundary conditions . - In flow on surfaces of the first kind, the 
flow picture is as shown in figure 9. Arhltrarily assigning a value 
tJtj on the suction sxirface, the value ijfj-j- on the pressure siirface of 

the next hlade is determined from the mass flow passing between them. 
These two values are used as the end values in eguatlon (l4l) for grid 
points next to the boundary. Outside the hlade region, however, the 
position of the dividing streamline is not known. Instead, there is 
the condition that the flow repeats itself or the ijr-value Increases 
hy when (p increases hy an amount equal to the pitch angle 

(2it divided hy number of blades). It is then convenient to draw any 
tjTO parallel lines up to the leading and trailing edges of the blade and 
consider only the grid points lying between the two reference lines. 

For the cp-derivative at a point for example, the required 

value is obtained from which is a pitch angle away from i{f^ 

(fig. 10), as 

ilfb = T|ff _ _ ilr^) (142) 

This relation is used between the inlet station 1-1 and the leading edge 
of the blade and between the exit station 2-2 and the trailing edge of 
the blade when the Sp surfaces are assumed to be surfaces of revolu- 
tion. For the general Sp s\irface where its deviation from the surface 
of revolution is considered, modification has to be made in places such 
as shown in the exit portion of figure 9. Because of the twist of the 
flow siirface, the dividing line from station 1-1 to the leading edge of 
the blade becomes two separate lines from the trailing edge of the blade 
to the exit station 2-2, accompanied by trailing vortices. Although the 
flow still repeats itself circumferentially every pitch angle, the use 
of equation (142) for the derivative at a point close to these lines 
will give inaccurate results. In these cases, it is better to use the 
end-point differentiation formulas to evaluate the derivatives. 
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At the station 1-1 sufficiently far upstream of the blades, the 
flow condition can he taken as uniform and the flow angle, equal to the 
given inlet angle. For the point h, the ijr value at point i 
upstream can he obtained from the given flow angle as follows: 


ijfi = TjrS + 


Thus, the required ijr value u; 
the values on that station, and only the ijr values downstream of 
station 1-1 will he Involved in the finite -difference expression (l4l) . 

An alternative method to take account of the inlet condition is 
as follows; If the first station 1-1 is chosen sufficiently far from 
the blades, the variation of the stream function upstream of the 
station 1-1 is linear in the circumferential direction. The value of 
the stream fimction, however, depends on the inlet angle. If solutions 
for a range of inlet angle are desired, they can he obtained by speci- 
fying a nimiber of sets of lineanly varying stream functions upstream of 
station 1-1 as fixed boundary values. !The slope of the streamlines 
obtained in the solution at the inlet then gives the value of the inlet 
angle. If, however, the solution for a certain specific inlet angle 
is desired, the streamline obtained in the solution must be adjusted 
according to that inlet angle, for example, as jk in figure 10 is 
adjusted to position gk, thereby obtaining an improved set of boundary 
values of the stream functions to be used in the next calculation. 

This method is, of course, not so accurate and convenient as the previ- 
ous method for obtaining a solution for a given inlet angle, but is 
desirable in the matrix solution because the inlet angle is then not 
involved in the matrix factorization, thereby making the same matrix 
factors usable for a range of inlet angle and Mach number. 

At the exit station far downstream of the blade, the same methods 
can be applied. For a blade having a sharp trailing edge, the Eutta- 
Joukowski condibion can be used and the correct exit angle far down- 
stream is the one that gives the flow at the trailing edge satisfying 
that condition. For round trailing edges, either the position of the 
stagnation point is assumed or some available enrplrlcal rule for the 
exit angle is used. If the calculation is made to con 5 )are with cert ain 
experimental results, the measured exit angle may be used. 


Sj, tan cxjj^ 




^ tan tan 0 ^) 


(143) 


pstream of station 1-1 can be replaced by 
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In flow along surfaces of the second kind, the houndary walls 
extend all the way to the inlet and exit stations with the T|r values 
given on the wal. 1 s (fig. H) . Across the inlet and exit, stations, the 
flow is considered to he uniform and pai’allel to the wadis so that the 
required lif value outside the station can he obtained hy an equation 
similar to equation (l43) . For the inlet station where the axial 
velocity is radially uniform and there is no radial or tangential 
velocity, ijr varies as the square of the radius. For the exit station 
with a certain radial gradient in fluid state, the radial variation of 
i{r can he determined from the corresponding radial variation in axieil 
velocity and density. 

Solution of finite-difference equations . - With the grid system 
and the degree of polynomial representation chosen and the houndary 
conditions taken into account, the prohlem remaining is the solution 
of the set of linear algebraic eq-uations (l4l) written for all interior 
grid points. For a small number of solutions with a given blade, the 
best method is the relaxation method (references 25, 33, 45, and 36). 

A modification of this method involving the use of higher-order differ- 
ences is suggested hy Fox (reference 46) . Formulas and tables of 
coefficients obtained in reference 26 enable the direct use of higher- 
degree polyncanials for problems with curved boundaries (reference 29). 
For the present flow problems, it is necessary to include a large 
domain to get to the boundary conditions that are given at places far 
from the blades, and the rise of higher-degree polynomials whenever it 
is applicable greatly reduces the numerical work. 

If a ntcmber of cases are to be solved for a given geometry (same 
blades for Sp surface and same hub and casing shapes for Sg sur- 
face), it is advantageous to solve the problem on a large-scale digital 
computing machine. If a high-speed digital machine is available, the 
simultaneous equations may be solved by Liebmann's Iterative process, 
which is the most simple to set up. For quicker results or when only a 
relatively slow-speed machine is available, the matrix process discussed 
in reference 26 is most suitable. In a calculation of the S2,m sur- 
face for a gas turbine and in a calculation of the Sp siorface of 
revolution for a centrifugal compressor, the coefficient matrices 
(about 400 and 200 interior grid points for the two problems, respec- 
tively, and the fourth-degree differentiation formula are used) were 
factorized into the lower and upper triangular matrices on an IBM CIEC 
and an IBM 604, respectively, in about 60 hours. The determination of 
i]r for a given set of values of K took only 2 hours on the CFEC for 
the gas-tvirbine problem. The gas-turbine problem was also worked out 
on an Univacj the factorization took only 11 minutes and the determina- 
tion of Tjr, 2.5 miinutes. The increasing availability of these high- 
speed large-scale digital calculating machines ^7111 render the suggested, 
method of solving the three-dimensional-flow problem a practical one. 
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General table for eval-uation of density from i^-derlvatives . - 
After the i{f values are obtained at the end of each cycle of calcula- 
tion, the velocity ccmponents are evaluated from the derivatives of ijr 
with respect to the coordinates, after the density is obtained as 
follows: From equations (46), (57), (74), (80), (85), (l05), and (113), 
the relation between h or p and ijr-derivatlves can be put into a 
ccamaon form as 



The quantities represented by X, as well as by e, t], and ^ for dif- 
ferent cases, ane given in the following table: 


Case 

Sunface 

Cooi^ilnates 

k 

X 

■ 

Equation 

■■ 


1 

(general) 

cp 

z 

1 

— w ^ 

2 r 

B 

(46) 

2 

Si (surface of 

cp 

z 

(1 + X2) 

0 

H 

(57) 


revolution) 







3 

Sq (cylindrical 


z 

1 

0 

1 



surface) 







4 

Si (general) 


r 

1 


I 

(74) 

5 

Si (surface of 


r 


0 

B 

(80) 


revolution) 



\ 


■ 


6 

Si (radial 


r 

1 

0 

B 

(85) 


plane) 





■ 


7 

S2 (Vui* given) 

r 

z 

1 

— W 2 

B 

(105) 

8 

32 given) 

r 

z 

1 

1 + 

0 

■ 

(113) 


With the Tir-derivatives evaluated, if an exact determination of h 
or p from the preceding equa,tion considering the vaxiatlon of specific 
heat with temperatmre is desired, the Keenan, and Kay gas tables (refer- 
ence 47) can be used. With two or three readings of h and p (or its 
reciprocal, specific volume), the correct value of h or p satisfying 
equation (144) is found. For most cases where the ten 5 >eratTxre range 
involved is not too large the use of an appropriate average value £f r, 
y, may give accurate enough resvilts. With the use of an average y, 
the density at any point in the flow field can be related to the inlet 
total value by equation (l2a) as 
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P 



s$^i-s* 




( 145 ) 


In order to mke out a general tatle for the calculation of density from 
the ijr-derivatives, the preceding equation is rewritten as 

2 


2 

2 



( 146 ) 



The functional relations between 2 sind $ are given in table I for y 
equal to 1.4 and 4/3, respectively. From the given inlet condition and 
the given X values, the variation of 

_ Lti 


% 




is 'first cc^^ted and plotted as an auxiliary graph or table as a func- 


tion of 


/ 2^2 \ 

^ ^ similar auxiliary graph or table is prepared 


for the variation of 


as a function of 

,. 2-_2 


I + I CD^r^ - X 


r-1 


(=|; ^ A 


H.- 


Anytime during the calculation, from the 


/ \ / 

value of 2 X 1 at each point (in general, X changes during 

successive improvements between S-j and So surfaces). 
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-2 


' I + J 


is read from the first graph or table and is combined with 


k 



+ 


1 

e aU 


and the entropy factor to obtain . The value of 2 is then read from 
tables I or II. After the value of 


I + 


2 2 
03 r 


- X 


% 


1 

f-1 


is read from the second curve or table, the density ratio is obtained. 


Hyperbolic Case 

In the hyperbolic case, the main problem is the solution of the 
following principal eq^uation, written in a common form for the two kinds 
of flow surface: 


K 

3^2 b 


+ 






+ M 




3ijr 


N at _ 
^ V an 


0 


(147) 


with the initial condition that and its normal derivative are given 
on a curve which is not a characteristic curve. From equation (147), 
the equation of the characteristic curve is 


The slopes of the characteristic curves are 


(148) 


A 




(149a) 
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^2 = = 

The coefficients J, K, L, and v and the independent variables 
Tj and ^ for the eight cases considered are given in the table on the 
following page. Using these values of and hj Aj_ and A 2 

also expressed in terms of the velocity cong)onents . Except for cases 2^ 
5, and 6, they can also be expressed in the usual trigonometric form, 
tan (X± p) . The values of X and p are also given. 

Changes of tlr-derivatives along characteristic curve . - inien the 
point nn the T]^ -plane moves along the image of the charac- 
teristic curve in the T]^-plane conre spending to a smal l change in 

d^, the change in t] is dr) = ^ d^. Because of these two s mal l 

changes, the change of any quantity q on the surface is (fig. 12) 


M + A ^ 
d^ 3C V ail 


(150) 


(151) 


Hence adong 


From equations 


d ^ a 3-^ Ai 

d^ 3^ “ 3^ at V 


d 3 ijr _ a 3 }Jf ^ ^ 
df an at 3 T 1 v 

(152) and (153), 


3^ _ 

a?5t “ At 


A §i - 4 

Ai 

(152) 

3t1 3^ 

a^2 

V 

3 3t _ a^ilr 

+ 

(153) 


a^aTi 

^ aTi^ 

3^ Ai ^ 
9T " V 9^2 


(154) 


_ A ^ A §1 +(^1 ^ 

3^2 dC a^ " V d^ 3 ti \v/ 3^2 


( 155 ) 
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Substituting equations (l54) and (l55) into equation (l4?) gives 




1 _d_ 8^ 
T? d^ St] 


+ (JAi^ 


2KAi + L) 


1 

~2~~2+^ 

3ri^ 


8tJt W 3ilf 

at V aq 


0 


(156) 


By virtue of equations (l49a) and (l49b), equation (l56) tecomeB 


■ 43i + ^§. + + J-!i = 0 

d^ a^ u d^ ail j a^ jv ari 

Similarly, along the second chsiracteristlc curve A 2 , 


(157a) 


d33, Al_da£ M23r H_|)t,o 

df u d^ ari J JV Sif 


(I57h) 


Starting from two points a and h a short distance apart on the 
initial curve, equations (l49a) and (l49h) give the tangent to the 
characteristic curves at these two points and equations (l57a) and 
(l57h) give the new value of and ^f/br] at the point of inter- 

section C of the two tangent lines (fig. ll) . The auxiliary equa- 
tions corresponding to the particular problem are then used to deter- 
mine other pertinent quantities at the point C. This process is to 
he carried step-hy-step downstream. 

Changes of fluid velocity and direction along characteristic 
curve . - When the characteristic cxirve hits the boundary wall, it is 
more convenient to express equations (l57a) and (l57b) in terms of the 
magnitude of the fluid velocity and the flow direction. In order to 
do this, the definitions of ■^-derivatives are first put in a common 
form for all cases as 


^ = bp Wf r 
9T1 ^ & 


_ 


bp Wt, € 


(158) 

(159) 


where £ equals 1 and r for the and Sg stirfaces, respectively. 

By the use of equation (45), 
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= - "bp Wt) € 


d In b 1 dh ■ ds * 

rsr~ ^2 af - ar 


d In ¥ti 


+ 


d la e 
d^ 


•) 


(160) 


d ^ 

dTa^ 


= Tap 




r 


d la b 1 dh ds * 

“45“"" ^ar - ^ + 


la W(; d In r\ 


(161) 


Sxibstitatlag equations 


(160) and (l6l) into equation (156) yields 





Ija b ^ 1 dh ds *\ (a 
^ ^ ;2 ^ 





+ 


A2 


d In r 

d5 




d la e . „ IT 

^r~ "t j - 




(162) 


Let 

¥^ = w sia X 
¥^ = -w cos X 

sad 


(163) 


h = I + 


CD^r^ 


- 1 («t‘ 


+ „ 2) 


= I + 


(jo^r^ 


I (v2 + ¥^2) 


(164) 


where ¥^ is equal to 0, ¥^, ¥^, 0, ¥^, aad ¥^ for cases 1 

to 8, respectively. By the use of equations (163), (l64), aad (l44), 
equation (162) can be written 


1/w^ -,\dw cosx + As slaX dv. d. L W*2\ 

^ \a2 / d^ Ag cos X - sin X d^ " a2 d^ V 2 " 2 / ~ 


^ ^ ^ i d la e , ydlnrN ds * 

A 2 cos X - sia X d^ “ ^ cosX + 

1 M siaX - ¥ cos X 
J A 2 COSX - sin X ~ ° 


(165) 
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A similar ejspression can "be obtained for the change in w and X 
along Ag by replacing Ag by in the preceding eqmtion. For 

cases 1, 3, 4, 6, and 7, A can be written as tan (X ±. where p is 

equal to sin"-^ — , throu^ purely trigonometric transformations as 
follows: 


A K +a/k2ZjL 
- J 


¥_ Wf ± a ^w2-a2 



_ sin X cos X ± sitL p cos p _ sin 2x ±. sin 2p 

~ cos 2X + cos 2p 


2 2 
cos X - sin p 


_ sin (X+p) cos (X-p) cos (X+p) sin (X-d) 

~ cos (x+p) cos (x-p) cos (x+p) cos (x-p) 


= tan (x+p) or tan (X-p) 

For these cases, equation (165) can then be written (compare refer- 
ence 30) : 


1 dw — . 

— + tan 

w d^ 


ax 


+ tan^ p 


a' 


1 d /_ ch^r^ 

-2 ar — 


^x 

2 


d In b . ds* . 1 „ d In e 

d^ d^ A cos X - sin X T ^ d^ 


A cos X 


d In r 




1 M sinX - W cos X _ q 
J A cos X - sin X L ~ 


(166) 


where the minus and plus signs on the second teim and subscripts 2 
and 1 for X in the last t^ro terms are used along characteristics Ai 

and Ag, respectively. Equations (l65) and {l66) are most useful when 

the characteristic hits the boundary wall. For a direct problem, the 
slope there is known from the given blade shape and for an Inverse or 
design problem, either the desired turning at the boundary or the 
velocity on the boundary is prescribed. With either dX or dw known, 
dw or dX is evaluated from eqiiation (162) or (l65) (only one charac- 
teristic eqxmtion is used at the wall) . For convenience of setup in 
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calciilation, this system can also he used for interior points. Except 
that more terms are involved in the present problem and that w takes 
different meanings in different cases, the procedure of calculation is 
very much the same as ordinary hvro-dlmensional flow described in refer- 
ences 32 and 30. 


CCNCLUDING ESMARKS 

A general theory of steady three-dimensional flow of a nonviscous 
fluid in subsonic and supersonic txirbomachines having arbitrary hub and 
casing shapes and a finite number of thick blades is presented. The 
solution of the three-dimensional direct and inverse problem is 
obtained by investigating a combination of flows on relative stream 
surfaces whose Intersection with a z-plane either upstream of or some- 
where inside the blade row form a circular arc or a radial line. The 
equations obtained to describe the fluid flow on these stream surfaces 
show clearly the several approximations involved in ordinary two- 
dimensional treatments. They also lead to a solution of the three- 
dimensional problem in a mathematically two-dimensional manner throu^ 
an iterative process. The equation of continuity is combined iflth the 
equation of motion in either the tangential or the radial direction 
through the use of a stream function defined on the surface, and the 
resulting equation is chosen as the principal equation for such flows. 
The character of this equation depends on the relative magnitude of 
the local velocity of soiind and a certain combination of velocity com- 
ponents of the fluid. A general method to solve this equation by both 
hand and machine computations when the equation is elliptic or hyper- 
bolic is described. The theory is applicable to both irrotational and 
rotational absolute flow at the inlet of the blade row and to both 
design and off -design operations. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, July 13, 1951 
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1 


1.0005007 

1.0010027 

1.0015061 

1.0020109 

1.0025171 


1.0030246 

1.0035335 

1.0040438 

1.0045555 

1.0050687 


1.0055833 

1.0060993 

1.0066168 

1.0071357 

1.0076561 


1.0081779 

1.0087012 

1.0092260 

1.0097522 

1.0102800 


1.0108093 

1.0113401 

1.0118724 

1.0124063 

1.0129417 


1.0134786 

1.0140171 

1.0145572 

1.0150989 

1.0156421 


1.0161669 

1.0167333 

1.0172814 

1.0178311 

1.0183825 


1.0217252 

1.0222882 

1.0228530 

1.0234194 

1.0239876 


1.0245575 

1.0251292 

1.0257027 

1.0262780 

1.0266550 


1.0274338 

1.0280145 

1.0285970 

1.0291813 

1.0297674 


0.99949955 

.99899830 

.99849616 

.99799314 

.99748922 


0.99698452 

.99647894 

.99597249 

.99546516 

.99495686 


0.99444770 

.99393768 

.99342669 

.99291486 

.99240207 


0.99188844 

.99137386 

.99085834 

.99034199 

.98962460 


0.98930629 

.98878706 

.98826690 

.98774573 

.98722365 


0.98670066 

.98617666 

.98565167 

.98512569 

.98459881 


0.98407094 

.98354210 

.98301217 

.98246128 

.98194932 


0.98141649 

.98088260 

.98054775 

.97981176 

.97927482 


0.97873675 

.97819773 

.97765759 

.97711652 

.97657433 


0.97603112 

.97548680 

.97494137 

.97439485 

.97384733 


0.97329872 

.97274893 

.97219805 

.97164610 

.97109308 


1 

I 


1.0303554 

1.0309453 

1.0315371 

1.0321307 

1.0327263 


1.0333238 
1.0339232 
1.0345246 
1.0351279 
1 .0357332 


1.0363404 

1.0369496 

1.0375609 

1.0381742 

1.0387895 


1.0394069 

1.0400263 

1.0406478 

1.0412714 

1.0418971 


1.0425249 

1.0431548 

1.0437869 

1.0444212 

1.0450576 


1.0456962 

1.0463370 

1.0469800 

1.0476253 

1.0482729 


1.0489227 

1.0495748 

1.0502292 

1.0508859 

1.0515449 


0.97053890 

.96998357 

.96942708 

.96886954 

.96831077 


0.96775086 

.96718983 

.96662757 

.96606419 

.96549961 


.96493392 

.96436702 

.96379885 

.96322948 

.96265894 


0.96208713 

.96151415 

.96093991 

.96036442 

.95978768 


0.95920970 

.95863049 

.95804996 

.95746812 

.95688506 


0.091 

.092 

.093 

.094 

.095 

1.0522063 

1.0528701 

1.0535362 

1.0542048 

1,0548758 

0.096 

.097 

.098 

.099 

.100 

1.0555492 

1.0562251 

1.0569034 

1.0575842 

1.0582675 

0.101 

.102 

.103 

.104 

.105 

1.0589534 

1.0596418 

1.0603327 

1.0610263 

1.0617225 

0.106 

.107 

.108 

.109 

.110 

1.0624213 

1.0631227 

1.0638268 

1.0645336 

1.0652431 


1.0659554 

1.0666704 

1.0673882 

1.0681087 

1.0688320 


1.0695582 

1.0702883 

1.0710192 

1.0717540 

1.0724918 


1.0732326 

1.0739763 

1.0747230 

1.0754728 

1.0762256 


1.0769815 

1.0777405 

1.0785026 

1.0792678 

1.0800362 


1.0808079 

1.0815828 

1.0823610 

1.0831424 

1.0839272 


1.0847155 

1.0855068 

1.0863017 

1.0671000 

1.0879018 


1.0887071 

1.0895160 

1.0903284 

1.0911444 

1.0919640 


1.0927673 

1.0936143 

1^0944450 

1.0952795 

1.0961178 


1.0969599 

1.0978058 

1.0986557 

1.0995095 

1.1003674 


1.1012292 

1.1020950 

1.1029650 

1.1038391 

1.1047174 


1.1055999 

1.1064866 

1.1073777 

1.1082731 

1.1091730 


0.93812555 

.93749672 

.93686627 

.93623430 

.93560073 


0.93496548 

.93432857 

.93369008 

.93304993 

.93240806 


0.93176447 

.93111924 

.95047232 

.92982361 

.92917321 


0.92852106 

.92786714 

.92721149 

.92655410 

.92589489 


0.92523380 

.92457092 

.92390616 

.92323964 

.92257118 


.92190089 

.92122868 

.92055458 

.91987858 

.91920061 


0.91852069 

.91783875 

.91715487 

.91646898 

.91578111 


0.91509116 

.91439916 

.91370512 

.91300896 

.91231070 


0.91161035 

.91090792 

.91020326 

.90949646 

.90878737 


0.90807618 

.90736280 

.90664708 

.90592913 

.90520888 


0.90448633 

.90376151 

.90303426 

.90230468 

.90157261 
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IEABI£ I - (SZCSRAL IWBUI £AB1£ - Contlim*d 
(a) X * - Concluded 



5 <p 

1 

X 

Z 


5 q > 

1 

z 

Z 


5 cp 

1 

T 

Z 

0.166 

1.1100772 

0.90083825 


0.221 

1.1683271 

0.85592468 


0.276 

1.2556050 

0.79642881 

.167 

1.1109859 

.90010143 


.222 

1.1695813 

.85500683 


.277 

1.2576894 

.79510887 

.168 

1.1118991 

.89936218 


.223 

1.1708444 

.85408445 


,278 

1.2598004 

.iBznesi 

.169 

1.1128170 

.89862035 


.224 

1.1721164 

.85315759 


.279 

1.2619389 

.79243139 

.170 

1.1137394 

.89787611 


.225 

1.1733975 

.85222612 


.280 

1.2641055 

.79107321 

0.171 

1.1146665 

0.89712932 


0,226 

1.1746880 

0.85128987 


0.281 

1.266301’2 

0.78970153 

.172 

1.1155983 

.89638000 


,227 

1.1759880 

.85034881 


.282 

1.2685266 

.78831602 

.173 

1.1165348 

.89562815 


.228 

1.1772974 

.84940305 


.283 

1.2707833 

.78691623 

.m 

1.1174762 

.89487364 


,229 

1.1786165 

.84845240 


.284 

1.2730715 

.78550184 

.175 

1.1184225 

.89411649 


.230 

1.1799455 

.84749677 


.285 

1.2753927 

.78407223 

0.176 

1.1193737 

0.89335670 


0,231 

1.1812846 

0.84653605 


0.286 

1.2777476 

0.78262718 

.177 

1.1203298 

.89259431 


.232 

1.1826339 

.84557021 


.287 

1,2801376 

.78116602 

.178 

1.1212910 

.89182915 


,233 

1.1839935 

.84459923 


.288 

1.2825638 

.77968831 

.179 

1.1222573 

.89106126 


,234 

1.1853636 

.84362300 


.289 

1.2850276 

.77819340 

.180 

1.1232287 

.89029064 


,235 

1.1867445 

.84264136 


.290 

1.2875299 

.77668099 

0.181 

1.1242053 

0.88951724 


0.236 

1.1881362 

0.84165435 


0.291 

1.2900725 

0.77515023 

.182 

1.1251871 

.88874108 


.237 

1.1895390 

.84066180 


.292 

1.2926566 

.77360066 

.183 

1.1261742 

.88796209 


.238 

1.1909530 

.83966370 


.293 

1.2952840 

.77203146 

.184 

1.1271667 

.88718022 


.239 

1.1923785 

.83865987 


.294 

1.2979563 

.77044196 

.185 

1.128 i 64,7 

.88639540 


.240 

1.1938156 

.83765030 


.295 

1.3006755 

.76883127 

0.186 

1.1291681 

0.88560773 


0,241 

1.1952645 

0.83663490 


0.296 

1.3034429 

0.76719893 

.187 

1.1301771 

.88481708 


.242 

1.1967255 

.83561351 


.297 

1.3062609 

.76554385 

.188 

1.1311917 

.88402346 


.243 

1.1981988 

.83458605 


.298 

1.3091314 

.76386526 

.189 

1.1322120 

.88322682 


.244 

1.1996846 

.83355242 


.299 

1.3120571 

.76216195 

.190 

1.1332381 

.88242709 


.245 

1.2011831 

.83251255 


.300 

1.3150400 

.76043314 

0.191 

1.1342699 

0.88162438 


0.246 

1.2026946 

0.83146628 


0.301 

1.3180832 

0.75867745 

.192 

1.1353077 

.88081848 


.247 

1.2042193 

.83041353 


.302 

1,3211890 

.75689398 

.193 

1.1363514 

.88000948 


.248 

1.2057574 

.82935423 


.303 

1.3243612 

.75508102 

.194 

1.1374012 

.87919724 


.249 

1.2073092 

.82828823 


. .304 

1.3276021 

V .75323774 

.195 

1.1384571 

.87838180 


.250 

1.2088750 

.82721539 


.305 

1.3309163 

.75136205 

0.196 

1.1395191 

0.87756318 


0.251 

1.2104551 

0.82613556 


0.306 

1.3343067 

0.74945288 

.197 

1.1405873 

.87674151 


.252 

1.2120497 

.82504868 


.507 

1.3377787 

.74750779 

.198 

1.1416619 

.87591607 


.253 

1.2136591 

.82395460 


.308 

1.3415356 

.74552558 

.199 

1.1427429 

.87508748 


.254 

1.2152836 

.82285320 


.309 

1.5449840 

.74350327 

.200 

1.1438304 

.87425548 


,255 

1.2169235 

.82174434 


.310 

li 3487277 

.74143951 

0.201 

1.1449245 

0.87342004 


0.256 

1.2185792 

0.82062783 


0.311 

1,3525748 

0.73933065 

.202 

1.1460252 

.87258116 


.257 

1.2202510 

.81950353 


.312 

1.3565299 

.73717505 

.203 

1.1471327 

.87173873 


.258 

1.2219392 

.81857132 


.313 

1.3606032 

.73496814 

.204 

1.1482470 

.87089276 


.259 

1.2236442 

.81723102 


.314 

1.3648001 

.73270804 

.205 

1.1493681 

.87004329 


.260 

1.2253663 

.81608251 


.315 

1.3691338 

.73038880 

0.206 

1.1504963 

0.86919011 


0.261 

1.2271059 

0.81492559 


0.316 

1,3736105 

0.72800841 

.207 

1.1516316 

.86833324 


.262 

1.2288634 

.81376010 


.317 

1.3782479 

.72555888 

.208 

1.1527741 

.86747265 


.263 

1.2306393 

.81258578 


.318 

1.3830531 

.72303804 

.209 

1.1539238 

.86660835 


.264 

1.2324338 

.81140261 


.319 

1.3880508 

.72043473 

.210 

1.1550810 

.86574015 


.265 

1.2342475 

.81021027 


.320 

1.3932489 

.71774684 

0.211 

1.1562457 

0.86486808 


0.266 

1.2360808 

0.80900860 


0.321 

1.3986833 

0,71495813 

.212 

1.1574179 

.86399217 


.267 

1.2379342 

.80779738 


.322 

1.4043617 

.71206727 

.213 

1.1585978 

.86311229 


.268 

1.2398081 

.80657644 


.323 

1.4103403 

.70904873 

.214 

1.1597855 

.86222840 


.269 

1.2417032 

.80534543 


.324 

1.4166227 

.70590426 

.215 

1.1609811 

.86134046 


.270 

1.2436198 

.80410428 


.325 

1.4233067 

.70258926 

0.216 

1.1621847 

0.86044843 


0.271 

1.2455586 

0.80285263 


0.326 

1.4303770 

0,69911639 

.217 

1.1633964 

.85955226 


,272 

1.2475200 

.80159036 


.327 

1.4380371 

.69539235 

.218 

1.1646164 

.85865183 


.273 

1.2495048 

.80031705 


.328 

1.4461815 

.69147614 

.219 

1.1658447 

.85774718 


.274 

1.2515134 

.79903260 


.329 

1.4554567 

.68706956 

.220 

1.1670816 

.85683812 


.275 

1.2535466 

.78773660 


.330 

1.4650468 

.68257205 
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TffllE I - (EIEHAL DEBSITT IfiBIS - Contimied 


6tp 

1 

Y 

£ 

0.001 

.002 

.003 

.004 

.005 

1.0005006 

1.0010026 

1.0015060 

1.0020107 

1.0025168 

0.99949965 

.99899840 

.99849626 

.99799333 

.99748952 

0.006 

.007 

.008 

.009 

.OlO 

1.0030242 

1.0035330 

1.0040432 

1.0045548 

1.0050678 

0.99698492 

.99647944 

.99597308 

.99546585 

.99495775 

0.011 

.012 

.013 

.014 

.015 

1.0055822 

1.0060980 

1.0066152 

1.0071339 

1.0076541 

0.99444879 

.99393896 

.99342827 

.99291663 

.99240404 

0.016 

.017 

.018 

.019 

.020 

1.0081756 

1.0086985 

1.0092230 

1.0097490 

1.0102764 

0.99189070 

.99137651 

.99086129 

.99034513 

.98982813 

0.021 

.022 

.023 

.024 

.025 

1.0108053 

1.0U3357 

1.0118676 

1.0124010 

1.0129359 

0.98931021 

.98879136 

.98827159 

.98775090 

.98722930 

0.026 

.027 

.028 

.029 

.030 

1.0134723 

1.0140104 

1.0145499 

1.0150910 

1.0156336 

0.98670679 

.98618318 

.98565876 

.98513335 

.98460705 

0.031 

.032 

.033 

.034 

.035 

1.0161778 

1.0167236 

1.0172710 

1.0178200 

1.0183707 

.98407975 

.98355148 

.98302222 

.98249199 

.98196069 

0.036 

.037 

.038 

.039 

.040 

1.0189229 

1.0194767 

1.0200322 

1.0205894 

1.0211482 

0.98142853 

.98089539 

.98036121 

.97982597 

.97928978 

0.041 

.042 

.043 

.044 

.045 

1.0217087 

1.0222708 

1.0228346 

1.0234001 

1.0239673 

0.97875255 

.97821438 

.97767518 

.97713494 

.97659369 

0.046 

.047 

.048 

.049 

.050 

1.0245363 

1.0251070 

1.0256794 

1.0262536 

1.0268295 

0.97605131 

.97550792 

.97496352 

.97441802 

.97387151 

0.051 

.052 

.053 

.054 

.055 

1.0274072 

1.0279867 

1.0285679 

1.0291510 

1.0297359 

0.97332392 

.97277523 

.97222556 

.97167471 

.97112279 


(b) r - 4/s 


z 



6<P 

1 

£ 

£ 

0.056 

.057 

.058 

.059 

.060 

1.0303226 

1.0309112 

1.0315016 

1.0320939 

1.0326880 

0.97056980 

.97001565 

.96946044 

.96890409 

.96834668 

0.061 

.062 

.063 

.064 

.065 

1.0332840 

1.0338819 

1.0344817 

1.0550835 

1.0356872 

0.96778814 

.96722846 

.96666766 

.96610563 

.96554249 

0.066 

.067 

.066 

.069 

.070 

1.0362928 

1.0369003 

1.0375099 

1.0381215 

1.0387350 

0.96497824 

.96441288 

.96384622 

.96327838 

.96270945 

0.071 

.072 

.073 

.074 

.075 

J.. 0393506 
1.0399681 
1.0405876 
1.0412093 
1.0418331 

0.96213924 

.96156796 

.96099550 

.96042169 

.95984664 

0.076 

.077 

.078 

.079 

.080 

1.0424589 

1.0430867 

1.0437167 

1.0443488 

1.0449830 

0.95927043 

.95869308 

.95811440 

.95753449 

.95695337 

■ 

1.0456194 

1.0462580 

1.0468987 

1.0475416 

1.0481868 

0.95637093 

.95578720 

.95520226 

.95461603 

.95402842 

0.086 

.087 

.088 

.089 

.090 

1.0488341 

1.0494837 

1.0501355 

1.0507896 

1.0514460 

0.95343963 

.95284948 

.95225807 

.95166530 

.95107119 

0.091 

.092 

.093 

.094 

.095 

1.0521047 

1.0527657 

1.0534290 

1.0540947 

1.0547628 

0.95047575 

.94987897 

.94928087 

.94868137 

.94808046 

0.096 

.097 

.098 

.099 

.100 

1.0554332 

1.0561061 

1.0567813 

1.0574590 

1.0581391 

0.94747825 

.94687456 

.94626958 

.94566314 

.94505533 

0.101 

.102 

.103 

.104 

.105 

1.0588217 

1.0595068 

1.0601944 

1.0608846 

1.0615773 

0.94444608 

.94383538 

.94322324 

.94260959 

.94199452 

0.106 

.107 

.108 

.109 

.110 

1.0622726 

1.0629704 

1.0636709 

1.0643740 

1.0650797 

0.94137795 

.94075997 

.94014041 

.93951938 

.93889687 



6cp 

1 

z 

£ 

O.IU 

.112 

.113 

.114 

.115 

1.0657881 

1.0664991 

1.0672128 

1.0679293 

1.0686486 

0.93827281 

.93764730 

.93702025 

.93639158 

.93576130 

0.116 

.117 

.118 

.119 

.120 

1.0693706 

1.0700954 

1.0708230 

1.0715534 

1.0722867 

0.93512951 

.93449612 

.93386115 

.93322461 

.93258641 

0.121 

.122 

.123 

.124 

.125 

1.0730229 

1.0737620 

1.0745040 

1.0752489 

1.0759968 

0.93194656 

.93130508 

.93066196 

.93001723 

.92937079 

0.126 

.127 

.128 

.129 

.130 

1.0767477 

1.0775017 

1.0782587 

1.0790187 

1.0797818 

0.92872267 

.92807278 

.92742122 

.92676800 

.92611304 

^3 1 • V 

Ip 1 S 3 

^1 1 • 9 
1 9 

1.0805481 

1.0813175 

1.0820901 

1,0828659 

1.0836449 

0.92545626 

.92479776 

.92413746 

.92347538 

.92281152 

0.136 

.137 

.138 

.139 

.140 

1.0844271 

1.0852126 

1.0860015 

1.0867937 

1.0875893 

0.92214590 

.92147843 

.92080904 

.92013783 

.91946473 

IWBl 

1,0883882 

1.0891906 

1.0899965 

1.0908058 

1.0916187 

0.91878982 

.91811295 

.91743414 

.91675347 

.91607079 

0.146 

.147 

.148 

.149 

.150 

1.0924351 

1.0932551 

1.0940787 

1.0949059 

1.0957368 

0.91538619 

.91469960 

.91401103 

.91332050 

.91262792 

0.151 

.152 

.153 

.154 

.155 

1.0965714 

1.0974098 

1.0982520 

1.0990979 

1.0999477 

0.91193332 

.91123662 

.91053784 

.90983706 

.90913413 

0.156 

.157 

.158 

.159 

.160 

1.1008014 

1.1016591 

1.1025207 

1.1033863 

1.1042559 

0.90842908 

.90772182 

.90701245 

.90630090 

.90558719 

0.161 

.162 

.163 

.164 

.165 

1.1051296 

1.1060074 

1.1068894 

1.1077756 

1.1086660 

0.90487125 

.90416308 

.90343263 

.90270990 

.90198491 
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TiffilE I - CEHERAI* DEHSITT !EABI£ - Conclxided 


6q> 

1 

1 

Z 

0.166 

.167 

.168 

.169 

.170 

1.1095607 

1.1104597 

1.1113630 

1.1122708 

1.1131830 

0.90125759 

.90052795 

.89979602 

.89906163 

.89832489 

0.171 

.172 

.173 

.174 

.175 

1.1140997 

1.1150210 

1.1159468 

1.1168773 

1.1178125 

0.89758574 

.89684410 

.89610006 

.89535350 

.89460442 

0.176 

.177 

.178 

.179 

.180 

1.1187524 

1.1196971 

1.1206466 

1.1216010 

1.1225603 

0.89385283 

.89309868 

.89234197 

.89158266 

.89082074 

0.181 

.182 

.183 

.184 

.185 

1.1235246 

1.1244940 

1.1254665 

1.1264481 

1.1274329 

0.89005617 

.88928887 

.88851887 

.88774618 

.88697075 

0.186 

.187 

.188 

.189 

.190 

1.1284230 

1.1294184 

1.1304192 

1.1314255 

1.1324372 

0.88619250 

.88541146 

.88462758 

.88384078 

.88305117 

0.191 

.192 

.193 

.194 

.195 

1.1334545 

1.1344774 

1.1355061 

1.1365405 

1.1375807 

0.88225862 

.88146313 

.88066458 

.87986306 

.87905851 

0.196 

.197 

.198 

.199 

.200 

1.1386269 

1.1396790 

1.1407372 

1.1418015 

1.1428719 

0.87825081 

.87744005 

.87662610 

.87580897 

.87498870 

0.201 

.202 

.203 

.204 

.205 

1.1439486 

1.1450317 

1.1461212 

1.1472172 

1.1483197 

0.87416515 

.87333827 

.87250807 

.87167452 

.87083762 

0.206 

.207 

.208 

.209 

.210 

1.1494289 

1.1505449 

1.1516677 

1.1527975 

1.1539342 

0.86999727 

.86915339 

.86830602 

.86745504 

.86660054 

0.211 

.212 

.213 

.214 

.215 

1.1550780 

1.1562290 

1.1573873 

1.1585530 

1.1597262 

0.86574240 

.86488057 

.86401501 

.86314567 

.86227249 

0.216 

.217 

.218 

.219 

.220 

1.1609069 

1.1620953 

1.1632916 

1.1644958 

1.1657080 

0.86139552 

.86051462 

.85962969 

.85874075 

.85784776 


(b) X n 4/5 - Concliided 
2 



6cp 

1 

z 

2 

0.221 

.222 

.223 

.224 

.225 

1.1669285 

1.1681568 

1.1693937 

1.1706391 

1.1718931 

0.85695068 

.85604946 

.85514399 

.85423424 

.85332015 

0.226 

.227 

.228 

.229 

.230 

1.1731558 

1.1744273 

1.1757079 

1.1769977 

1.1782966 

0.85240170 

.85147884 

.85055140 

.84961933 

.84868275 

0.231 

.232 

.233 

.234 

.235 

1.1796049 

1.1809228 

1.1822504 

1.1835879 

1.1849354 

084774148 

.84679540 

.84584450 

.84488866 

.84392786 

0.236 

.237 

.238 

.239 

.240 

1.1862930 
1.1876610 
1-. 1890394 
1.1904285 
1.1918285 

0.84296207 

.84199111 

.84101502 

.84003365 

.83904689 

0.241 

.242 

.243 

.244 

.245 

1.1832396 
1.1946617 . 
1.1960954 
1.1975406 
1.1989976 

0.83805472 

.83705705 

.83605371 

.83504476 

,.83403003 

0.246 

.247 

.248 

.249 

.250 

1.2004666 

1.2019479 

1.2054416 

1.2049480 

1.2064673 

0.83300943 

.83198282 

.83095017 

.82991133 

.82886623 

0.251 

.252 

.253 

.254 

.255 

1.2079997 

1.2095455 

1.2111050 

1.2126763 

1.2142658 

0.82781478 

.82675683 

.82569224 

.82462101 

.82354292 

0.256 

.257 

.258 

.259 

.260 

1.2158677 

1.2174844 

1.2191161 

1,2207632 

1.2224259 

0.82245790 

.82136576 

.82026642 

.81915969 

.81804549 

0.261 

.262 

.263 

.264 

.265 

1,2241045 

1.2257994 

1.2275110 

1.2292396 

1.2309855 

0.81692372 

.81579417 

.81465665 

.81351105 

.81255725 

0.266 

.267 

.268 

.269 

.270 

1.2327492 

1.2345310 

1.2363514 

1.2381508 

1.2399896 

0.8U19501 

.81002421 

.80884462 

.80765606 

.80645838 

0.271 

.272 

.273 

.274 

.275 

1.2418483 

1.2437273 

1.2456273 

1.2475485 

1,2494916 

0.80525133 

.80403478 

.80280835 

.80157204 

.80032551 


6tp 

1 

2 

2 

0.276 

.277 

.278 

.279 

.280 

1.2514571 

1.2534457 

1.2554579 

1.2574944 

1.2595556 

0.79906854 
.79780081 
. 79652213 
.79523217 
.79393081 

0.281 

.282 

.283 

.284 

.285 

1,2616424 

1.2637555 

1,2658957 

1.2680636 

1.2702602 

0.79261762 

.79129230 

.78995450 

.78860398 

.78724028 

0.286 

.287 

.288 

.289 

.290 

1.2724861 

1.2747424 

1.2770299 

1.2793498 

1.2817028 

0,78586320 

.78447222 

.78306702 

.78164705 

.78021207 

0.291 

.292 

,293 

.294 

.295 

1.2840902 

1.2865130 

1.2889726 

1.2914700 

1.2940068 

0.77876149 
.77729490 
.77581168 
.77431144 
. .77279347 

0.296 

.297 

.298 

.299 

.300 

1.2965843 

1.2992041 

1.3018675 

1.3045765 

1.3073325 

0.77125722 

.76970200 

.76812732 

.76653228 

.76491635 

0,301 

.302 

.303 

.304 

.305 

1.3101379 

1.3129943 

1.3159043 

1.3188697 

1.3218935 

0.76327843 

.76161793 

.75993368 

.75822502 

.75649059 

0.306 

.307 

.308 

.309 

.310 

1.3249779 
1 .^281263 
1.3313412 
1.3346266 
1.3379854 

0.75472957 

.75294044 

.75112225 

.74927324 

.74739231 

0.311 

.312 

.313 

.314 

.315 

1.3414224 

1.3449408 

1.3485465 

1.3522432 

1.3560384 

0.74547734 
. 74352715 
.74153913 
.73951195 
.73744224 

0.316 

.317 

.318 

.319 

.320 

1.3599361 

1.3639457 

1.3680720 

1.3723269 

1.3767161 

0.73532867 

.73316702 

.73095568 

.72868935 

.72636617 

0.321 

.322 

.323 

.324 

.325 

1.3812551 

1.3859506 

1.3908236 

1.3958818 

1.4011550 

0.72397923 

.72152644 

.71899844 

.71639304 

.71369691 

0.326 

.327 

.328 

.329 

.330 

1.4066513 

1.4124156 

1.4184546 

1.4248420 

1.4315749 

0.71090824 

.70800691 

.70499260 

,70183220 

.69853139 
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Figure 1. - Belative stream surface 
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Figure 3. - Intersecting and Sg surfaces In a tlade rovr 
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INLET GUIDE VANE ROTOR STATOR 



Figure 5. - Mean stream surfaces far inlet stage of axial- flow compressor. 







MCA TN 2604 


89 












NACA 03J 2604 


93 


T 1 or r 



Figure 12. - Characteristic system for hyperbolic -case. 


NACA - Langley Field, Va. 


